Los Alamos Electronic Archives: physics/990606^ 



0^ ■ IFUG 

^ ' La verdad os hard libres 

CURSO: MECANICA CLASICA 



O 



EDITOR: HARET C. ROSU 
rosu@ifug3.ugto.mx 



f:--ha 



(N 
> 

\o 
o 
\o 

O I curso de maestria 

3^ ■ (graduate course) 

O 

•i-H 



X 



& 



Copyright © 1999 H.C. Rosu 

Leon, Guanajuato, Mexico 

vl: Junio de 1999. 



INDICE DE CONTENIDO 

1. LOS PRINCIPIOS DE MINIMO ... 3. 

2. MOVIMIENTO BAJO FUERZAS CENTRALES ... 20. 

3. CUERPO RIGIDO ... 34. 

4. OSCILACIONES PEQUENAS ... 55. 

5. TRANSFORMACIONES CANONICAS ... 73. 

6. PARENTESIS DE POISSON ... 83. 

7. LAS ECUACIONES DE HAMILTON- JACOBI ... 86. 

8. VARIABLES ACCION-ANGULO ... 94. 

9. TEORIA CANONICA DE PERTURBACIONES ... 101. 

10. INVARIANTES ADIABATICOS ... 116. 

11. MECANICA DE SISTEMAS CONTINUOS ... 120. 

Estudiantes colaboradores: 

Eri Mena (1) 

Julio Lopez (2,3) 

Alberto Juarez (4) 

Mario Ranferi Gutierrez (5-8) 

Zaida Urrutia (9,10) 

Monica Beltran (11) 

Gran parte de la responsabilidad del idioma pertenece a los estudiantes. 



1. LOS PRINCIPIOS DE MINIMO 



Prologo: La historia de los principios de "miniino" en la fisica es larga 
e interesante. La investigacion de tales principios se argumenta sobre la 
idea de que la Naturaleza actiia siempre de tal forma que determinadas 
cantidades de importancia resultan ser siempre minimizadas cuando tiene 
lugar un proceso fisico. La base matematica de estos principios es el calculo 
variacional. 
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1. Introduccion 

La experiencia ha demostrado que, cuando sea posible despreciar los efec- 
tos relativistas, el moviniiento de una particula dentro de un sistema de 
referenda inercial queda correctaniente descrito mediante la ecuacion de 
Newton F = dp/dt. Cuando suceda que la particula no haya de ejecutar 
un movimiento complicado y se utilicen coordenadas rectangulares para de- 
scribirlo, generalmente las ecuaciones de movimiento seran relativamente 
sencillas; aliora bien si no se verifica ninguna de estas condiciones, las ecua- 
ciones pueden hacerse bastante complicadas y dificiles de manejar. 
Cuando una particula esta limitada a moverse sobre una superficie dada, 
deben existir ciertas fuerzas (llamadas fuerzas de ligadura) que mantengan 
a la particula en contacto con dicha superficie. Con el fin de facilitar al- 
gunos problemas de indole practico que aparecen al aplicar las formulas 
de Newton a ciertos problemas, pueden desarrollarse otros procedimientos. 
Esencialmente, todos estos procedimientos de abordar los problemas son a 
posteriori, puesto que sabemos de antemano que hemos de obtener resul- 
tados equivalentes a las formulas de Newton. Entonces, no es necesario 
formular una nueva teoria de la mecanica -la teoria de Newton es suficien- 
temente correcta- para efectuar una simplificacion, sino que basta con idear 
otro metodo que nos permita abordar problemas complicados de forma gen- 
eral. El principio de Hamilton contiene un metodo de este caracter y las 
ecuaciones de movimiento que resultan de la aplicacion del mismo se llaman 
ecuaciones de Lagrange. 

Si las ecuaciones de Lagrange han de constituir una descripcion adecuada 
de la dinamica de las particulas, deberan ser equivalentes a las ecuaciones 
que resulten de las formulas de Newton. Por otra parte, el principio de 
Hamilton es de aplicacion a una amplia gama de fenomenos fisicos (en es- 
pecial los relatives a campos) con los que generalmente no se relacionan las 
ecs. de Newton. Es seguro que cada una de las consecuencias que pueden 
extraerse del principio de Hamilton fue deducida primero, al igual que las 
ecs. de Newton, relacionando entre si hechos experimentales. El princi- 
pio de Hamilton no nos proporciona teoria fisica nueva alguna, pero nos ha 
permitido unificar satisfactoriamente muchas teorias separadas, partiendo 
de un postulado fundamental sencillo. Ello no constituye un ejercicio de 
habilidad, puesto que el objetivo de la fisica no es unicamente dar una for- 
mulacion matematica precisa para los fenomenos observados, sino tambien 
describir sus efectos con ahorro de postulados fundamentals y de la manera 
mas unificada posible. 



El primero de los principios de minimo se desarrollo en el campo de la optica 
por Heron de Alejandria hace casi dos mil afios. Encontro que la ley que rige 
la reflexion de la luz puede obtenerse admitiendo que un rayo luminoso, que 
viaje de un punto a otro reflejandose en un espejo piano, recorre siempre el 
camino mas corto. No obstante, el principio del camino mas corto de Heron 
no puede proporcionar una expresion correcta de la ley de la refraccion. En 
1657, Fermat formulo nuevamente el principio postulando que los rayos lu- 
minosos viajan siempre de un punto a otro de un medio siguiendo el camino 
que requiera el menor tiempo. Este principio del tiempo minimo de Fer- 
mat conduce inmediatamente, no solo a la ley correcta de la reflexion, sino 
tambien a la ley de la refraccion de Snell. 

Los estudios acerca de los principios de minimo continuaron y, en la ultima 
parte del siglo XVH, Newton, Leibniz y los hermanos Bernoulli iniciaron el 
desarrollo del calculo variacional. En afios posteriores este principio recibio 
de Lagrange una base matematica solida (1760). En 1828, Gauss desarrollo 
un metodo para estudiar la mecanica mediante su principio de ligadura 
minima. En sendos trabajos, publicados en 1834 y 1835, Hamilton expuso 
el principio dinamico de la minima accion sobre el cual es posible formular 
toda la mecanica y a decir verdad, la mayor parte de la fisica. 
Accion es una magnitud de dimensiones longitud por I'mpetu o bien energia 
por tiempo. 

2. Principio de minima accion 

La formulacion mas general de la ley del movimiento de los sistemas mecanicos 
es el principio de minima accion (a de Hamilton). Segiin este principio, todo 
sistema mecanico esta caracterizado por una funcion definida: 

L(qi,q2,...,qs,qi,q2,qs,t 



o mas brevemente L {q,Q,t\, y el movimiento del sistema satisface la sigu- 

iente condicion: Supongamos que en los instantes t = ti y t = ^2 el sistema 
ocupa posiciones dadas, caracterizadas por los dos conjuntos de valores de 
las coordenadas q^^' y q^'^'; el sistema se mueve entre estas posiciones de 
manera que la integral 

S = f\ [q,q,t^ dt (1) 

tome el menor valor posible. La funcion L se llama Lagrangiana del sistema, 
y la integral (1) la accion. La funcion de Lagrange no contiene mas que q 



y (7, y no las derivadas superiores, eso es debido al hecho que el estado 
mecanico de un sistema esta completamente definido por sus coordenadas y 
sus velocidades. 

Establezcamos ahora las ecuaciones diferenciales que determinan el mmimo 
de la integral (1). Por simplicidad empecenios suponiendo que el sistema no 
tiene mas que un solo grado de libertad, de manera que hace falta determinar 
una sola funcion q (t). Sea precisamente q = q{t) la funcion para la cual S 
es un mmimo. Esto significa que S crece cuando se sustituye q (t) por una 
funcion cualquiera 

q{t) + 6q{t), (2) 

donde 5q (t) es una funcion que es pequeiia en todo el intervalo de tia t2 [se 
le llama variacion de la funcion q (t)]. Puesto que para t = ti y t = t2 todas 
las funciones (2) deben tomar los mismos valores q^^' y q^'^' , se tiene: 

6q ih) = dq (t2) = 0. (3) 

Lo que varia S cuando se reemplaza q por q + 5q esta dado por: 

/ L(q + 5q,q+6q,tjdt- L(q,q,tjdt. 

El desarrollo en serie de esta diferencia en potencias de 5q y 5 Q comienza 
por terminos de primer orden. La condicion necesaria de mmimo (en general 
extremal) de S es que el conjunto de estos terminos se anule; Asi el principio 
de minima accion puede escribirse: 

SS = S f\(q,q, t] dt = 0, (4) 



>ti 
o, efectuando la variacion: 

/ T^Sq + — S q]dt = . 
Jt2 \oq Qq J 

Teniendo en cuenta que 6 q= d/dt {5q), integramos el sugundo termino por 
partes y se obtiene: 



5S 



dL ^ 

^6q 

dq . 



*' "*i fdL d dL 



"Uw-s^j^""*^"' <^' 



En virtud de las condiciones (3), el primer termino de esta expresion desa- 
parece. Queda una integral, la cual debe anularse para todo valor de 5q. 
Esto es solamente posible si el integrando es identicamente nulo, y conse- 
cuentemente se obtiene la ecuacion: 

dL d dL _ 
dq dt Qq 

Si hay varies grades de libertad, las s funciones diferentes qi{t) deben variar 
independienteniente. Es evidente que entonces obtenemos s ecuaciones de 
la forma: 

-j-^ = (. = 1.2,.....) (6) 

Estas son las ecuaciones diferenciales buscadas; en Mecanica se les llama 
ecuaciones de Lagrange. Si se conoce la lagrangiana de un sistema mecanico 
dado, entonces las ecuaciones (6) establecen la relacion entre las acelera- 
ciones, las velocidades y las coordenadas, es decir, son las ecuaciones del 
movimiento del sistema. Desde un punto de vista matematico, las ecua- 
ciones (6) forman un sistema de s ecuaciones diferenciales de segundo orden 
con s funciones desconocidas qi{t). La solucion general del sistema contiene 
2s constantes arbitrarias. Para determinarlas y, por lo tanto, para definir 
completamente el movimiento del sistema mecanico, es necesario conocer las 
condiciones iniciales que caractericen el estado del sistema en un instante 
dado, por ejemplo los valores iniciales de las coordenadas y de las veloci- 
dades. 

3. Principio de D'Alembert 
Desplazamiento virtual (infinitesimal) de un sistema es el cambio de con- 
figuracion de este a consecuencia de una variacion infinitesiamal arbitraria 
de las coordenadas 5ri, compatible con las fuerzas y ligaduras impuestas al 
sistema en el instante dado t. Se llama virtual al desplazamiento para dis- 
tinguirlo del desplazamiento real del sistema que tiene lugar en un intervalo 
de tiempo dt, durante el cual puede variar las fuerzas y ligaduras. 
Las ligaduras (o restricciones) introducen 2 tipos de dificultades en la solucion 
de problemas mecanicos: 

(1) No todas las coordenadas son independientes. 

(2) Las fuerzas de ligadura por lo general no se conocen a priori, son parte 
de las incognitas del problema y ban de obtenerse a partir de la solucion 
buscada. 



En el caso de ligaduras holonomicas la dificultad (1) se salva introduciendo 
coordenadas independientes {qi,q2,...,Qm, donde m es el numero de grades 
de libertad). Esto es, si hay m ecuaciones de ligadura y 3A^ coordenadas 
(xi, ...,X3n), podemos eliminar esas n ecs. introduciendo las variables inde- 
pendientes {qi,q2i ■■! iQn) niediante una transformacion de la forma 

xi = fi{qi,...,qm,t) 



X3N = fsNiqi, ■■■,qn,t) , 

donde n = 3N — m. 

Para librarnos del problenia (2) necesitamos formular la mecanica de modo 

que las fuerzas de ligadura NO APAREZCAN en la solucion de problemas. 

Esta labor constituye la esencia de lo que llamarenios "M Principio de Tra- 

bajos Virtuales" . 

Trabajo Virtual: Supongamos que un sistema de N particulas se describe 

por 3A^ coordenadas {xi,X2, ■■■,X3n) Y sean FiF2^...F^n las componentes 

de las fuerzas que actiian sobre cada uno. Si las particulas del sistema ex- 

perimentan desplazamientos infinitesimales e instantaneos 8x1^8x2, ...,8x31^ 

debido a dichas fuerzas, entonces el trabajo realizado es: 

3N 

8W = J2 Pj^xj . (7) 

i=i 

Dichos desplazamientos se llaman desplazamientos virtuales y 5W se llama 
trabajo virtual; (7) tambien puede escribirse como: 

N 

6W = J2'^»-Sr . (8) 

Q = l 

Fuerzas de ligadura o de restriccion: ademas las fuerzas externas F^ 
las particulas pueden estar sujetas a fuerzas de ligadura F^. 
Principio de trabajo virtual: Sea Fq, la fuerza que actua sobre la a- 
esima parti'cula, si separamos Fq, en la contribucion de origen externo F^ 
y la ligadura Rq, 

F, = Fi^) + R, . (9) 



Si el sistema esta en equilibrio entonces 

F, = F(f) + R, = . (10) 

Asi que el trabajo virtual debido a todas las posibles fuerzas F^ es: 

N N 

W=Y,Fa-Sr^ = J2 (^i'^ + Ra) • <5r„ = . (11) 

Q=l a=l 

Si el sistema es tal que sus fuerzas de ligadura no producen trabajo virtual 
entonces de (11) concluimos que: 

N 

5:F(f).5r„ = 0. (12) 

Ya hechan las definiciones anteriores podemos llegar a lo que es el principio 
de D'Alembert. La ecuacion del movimiento es segiin Newton: 

y puede escribirse de la forma 

Fa- Pa= , 

que dice que las particulas del sistema estaran en equilibrio bajo una fuerza 

igual a la real mas una fuerza invertida — Pj . En vez de (12) podemos 
escribir inmediatamente 

E(Fa-Pa -(Jr, =0 (13) 

y haciendo la misma descomposicion en fuerzas aplicadas y fuerzas de lig- 
adura (fa), resulta: 

N . s N 

a=l ^ ^ a=l 

Limitemonos de nuevo a sistemas para los cuales el trabajo virtual de las 
fuerzas de ligadura sea nulo y obtendremos 

E F(f)-pJ.<5ra = 0, (14) 



a=l 



que constituye el principio de D'Alembert. Ahora esta ecuacion anterior 
aiin no tiene forma litil para proporcionar las ecuaciones de movimiento del 
sistema, por lo que debemos transfomar el principio en una expresion que 
contenga desplazamientos virtuales de las coordenadas generalizadas, las 
cuales son entonces independientes entre si esto implica que se podran hacer 
cero los coeficientes de Jq^ y la velocidad en termino de las coordenadas 
generalizadas es: 

dt ^ dqk at 

Analogamente, el desplazamiento virtual arbitrario dr^ se puede relacionar 
con los desplazamientos virtuales 5<\j mediante 

^r« = E ^% • 

Entonces en funcion de las coordenadas generalizadas, el trabajo virtual de 
las Fq sera: 

donde las Qj son las llamadas componentes de la fuerza generalizada, las 
cuales se definen en la forma 

Ahora si vemos a la ec. (14) como: 

E P -^^c, = ^mara -dva (16) 

y sustituyendo en esta ultima ec. los resultados anteriores podemos ver que 
(16): 



^r^ \ d ( dVa \ dVa \ v;- 



dt[dqj dq,j ^' 

(17) 



6qj = 
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Las variables qj pueden ser un sistema cualquiera de coordenadas para de- 
scribir el movimiento del sistema. Sin embargo, si las ligaduras son holono- 
mas, sera posible encontrar sistemas de coordenadas qj independientes que 
contengan implicitamente las condiciones de ligadura en las ecuaciones de 
transformacion Xi = fi es que se anulen por separado los coeficientes: 

En total hay m ecuaciones. Las ecuaciones (18) suelen llamarseles ecuaciones 
de Lagrange, si bien esta designacion se reserva frecuentemente para la forma 
que toman estas ecuaciones, cuando las fuerzas se derivan de un potencial 
escalar V 

Fa = -ViV. 

Entonces Qj puede escribirse como: 

dqj 

Las ecuaciones (18) pueden escribirse tambien en la forma: 

dt\dqjj dqj 

y definiendo la funcion Lagrangiana L, en la forma L = T — V se obtiene 

d ( dL\ dL 



dt \dqjl dqj 
Estas son las Ecuaciones de Lagrange. 



. (20) 



4. - Espacio Fasico 

En la interpretacion geometrica de los fenomenos mecanicos se hace fre- 
cuentemente uso del concepto de espacio fasico, es un espacio de 2s di- 
mensiones cuyos ejes coordenados corresponden a los s coordenadas gener- 
alizadas y a los s impetus del sistema mecanico considerado. Cada punto 
en este espacio corresponde a un estado definido del sistema. Cuando el 
sistema se mueve, el punto representativo en el espacio fasico describe una 
linea denominada trayectoria fdsica. 
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5. - Espacio de Configuraciones 

El estado de un sistema compuesto de n particulas y sometido a m ligaduras 
que relacionen entre si a algunas de las 3n coordenadas rectangulares queda 
especificado por conipleto mediante s = 3n — m coordenadas generalizadas. 
Es posible, pues, representar el estado de tal sistema por un punto de un 
espacio de s dimensiones que llamamos espacio de configuraciones, corre- 
spondiendo cada una de las diemensiones de este espacio a una de las qj . La 
historia, o evolucion a traves del tiempo, del sistema, estara representada 
por una curva del espacio de confuguraciones, cada uno de cuyos puntos 
represaentara la configuracion del sistema en un instante determinado. 

6. - Ligaduras 

Es necesario tener en cuenta las ligaduras que limitan el movimiento del 
sistema. Las ligaduras pueden clasificarse de divesas maneras. En el caso 
general en que la ecuaciones de ligadura puedan expresarse como: 

/ , Cqj Qi — U , 

i 

donde las Cai son funciones de las coordenadas solamente (el indice a numera 
las ecuaciones de ligadura). Si los primeros miembros de estas ecuaciones 
no son derivadas totales con respecto al tiempo de funciones de las coor- 
denadas, estas ecuaciones no pueden ser integradas. En otras palabras, no 
pueden reducirse a relaciones entre las coordenadas solamente, que podrian 
utilizarse para expresar la posicion e los cuerpos por un niimero menor de 
coordenadas, correspondiente al mimero real de grados de libertad. Tales 
ligaduras se llaman no holonomicas (en oposicion a las ligaduras anteriores 
son las Uamadas holonomicas que relacionan solamente las coordenadas del 
sistema^. 

7. Ecuaciones de movimiento de Hamilton 

La formulacion de las leyes de la Mecanica con la ayuda de la Lagrangiana, 
presupone que el estado mecanico del sistema esta determinado dando sus 
coordenadas y velocidades generalizadas. Sin embargo, este no es el linico 
metodo posible; la descripcion del estado de un sistema en funcion de sus 
coordenadas e impetus generalizados presenta un cierto mimero de ventajas. 
El paso de un conjuto de variables independientes a otro puede realizarse 
mediante lo que se llama en matematicas tranformacion de Legendre. En 
este caso la transformacion toma la siguinte forma, donde la diferencial total 
de la Lagrangiana como funcion de las coordenadas y de las velocidades es: 

i ^^* i dqi 

12 



la cual puede escribirse como: 

dL = ^Pidqi+^Pidqi , (21) 

i i 

donde ya sabemos que las derivadas dL/d qi, son por definicion, los impetus 
generalizados y dL / dqt =pi por las ecuaciones de Lagrange. El segundo 
termino de la ec. (21) puede escribirse en la forma 

^Pid Qi=d(^piqA -^qi dqi . 



Llevando la diferencial total di^piqA al primer miembro, y cambiando 
los signos, se obtiene de (21): 

La cantidad bajo el signo de la diferencial es la energia del sistema expre- 
sada en funcion de las coordenadas y de los impetus y se llama funcion de 
Hamilton o Hamiltoniana del sistema: 

H{p,q,t)=Y,Pi%-L. (23) 

i 

Entonces de la ec. (22) 

dH = -Y,P, dqi + Y.Pi'^i 

en lo cual las variables independientes son las coordenadas y los impetus, se 
obtienen las ecuaciones 

■ dH ■ dH , ,, 

dpi dqi 

Estas son las ecuaciones de movimiento en las variables q y p y se llaman 
ecuaciones de Hamilton. 

8. Leyes de conservacion 
8.1 Energia 

Consideremos primero el teorema de conservacion que resulta de la homo- 
geneidad del tiempo. En virtud de esta homogeneidad, la Lagrangiana de un 
sistema cerrado no depende explicitamente del tiempo. Entonces la derivada 
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total respecto al tiempo de la Lagrangiana (no dependiente explicitamente 
del tiempo) puede escribirse como: 

dL .r-^ dL ■ .r-^ dL ■■ 
dt ^dq^ ^ Q q. 

y de acuerdo a las ecs. de Lagrange podemos reescribir la ec. anterior como: 




-L =0 . 

^q^ J 

De donde se deduce que la magnitud 



E = Y,q,^-L (25) 

i dQi 

permanece constantc durante el movimiento de un sistema cerrado, es decir 
es una integral del movimiento. A esta magnitud se le llama energia E del 
sistema. 
8.2 Impetu 

De la homogeneidad del espacio se deduce otro teorema de conservacion. 
En virtud de dicha homogeneidad, las propiedades mecanicas de un sistema 
cerrado no varian por un desplazamiento paralelo de todo el sistema en el es- 
pacio. Consideremos un desplazamiento infinitesimal e (i.e., los vectores de 
posicion r^ se convierten en r^ + e) y busquemos la condicion para que la la- 
grangiana no varie. La variacion de la funcion L, consecuencia de un cambio 
infinitesimal en las coordenadas (permaneciendo constantes las velocidades 
de las particulas), esta dado por: 

a " a " 

extendiendo la suma a todas las particulas del sistema. Como e es arbitrario, 
la condicion JL = es equivalente a 

\-^ dL , , 
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y en virtud de las ecs. de Lagrange ya mencionadas 

d f dL 



E 



dt \5Ve 



d y—^ dL 



dt 



. 



Llegamos asi a la conclusion de que en un sistema mecanico cerrado, la 
magnitud vectorial (llamada impetu) 



p-E 



dL 

dVa 



permanece constante durante el movimiento. 
8.3 Momento angular 6 cinetico 

Estudiemos ahora el teorema de conservacion que infiere de la isotropia del 
espacio. Consideremos una rotacion infinitesimal del sistema, y busquemos 
la condicion para que la Lagrangiana no varie. 

Llamaremos vector de rotacion infinitesimal (50 al vector cuyo modulo es 
igual al angulo de rotacion 6(p y cuya direccion coincide con el eje de rotacion. 
Consideremos primero el incremento en el vector de posicion correspondiente 
a una particula del sistema, tomando un origen de coordenadas situado en 
el eje de rotacion. El desplazamiento lineal extremo del vector de posicion 
en funcion de angulo es 



\6r\ 



r sm 



(ver fig.). La direccion del vector 6r es perpendicular al piano definido por 
r y 6(j), y por tanto, 



5r 



X r 



(27) 
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La rotacion del sistema no solamente modifica la direccion de los vectores 
de posicion sino tambien las velocidades de las particulas, transformandose 
en todos los vectores segiin la misnia ley. El increniento de velocidad con 
respecto a un sistema fijo de coordenadas sera; 

(5v = (50 X V . 

Llevemos estas expresiones a la condicion de que la Lagrangiana no varia 
por la rotacion: 

^ \dra dva J 

y sustituyendo, por definicion las derivadas dL/dva por p^ y las derivadas 
dL/dva de acuerdo con las ecs. de Lagrange, por P^; obtenemos 

J2 I Pa ^(px r^ + pa-6(l)XvA =0 , 

a ^ ^ 

o permutando circularmente los factores y sacando Scj) fuera del signo suma: 

puesto que Sep es arbitrario, resulta 

d 



dt 



^r„xp„ = 



y se concluye que en el movimiento de un sistema cerrado se conserva la 
magnitud vectorial (llamada momento angular 6 momento cinetico) 

M = ^r„xp„ • 



9.- Aplicaciones del principio de Accion 
a) Ecuaciones de movimiento 

Hallar las ecuaciones del movimiento para una masa pendular suspendida 
de un resorte, por aplicacion directa del principio de Hamilton 
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Para el pendulo de la figura la Lagrangiana es de la forma 



1.2 -2 i 
L = -m{r +r ) + mgr cos k{r — Vo)"^ , 



por lo tanto 



t2 ri2 

5Ldt = / 
ti Jti 



.2 



m [r S r +r +r 9 6 O] + mg5r cos — mgr69 sin 9 — k{r — ro)Sr 



dt 



mr6rdt = mr d{5r) = dimr 5r\ — m5r r dt 



Igualmente 



d I mr"^ 
mr^0^5 9 dt = d( mr^ 9 S 0] - ^0—^ ^-dt 



d mr^ 6 0] - 69 { mr^ +2mr r0 dt 



Por tanto la integral anterior se escribe como 
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t2 



m r —mr 



mg cos 6 + k {r — Tq) > + s fnr +2mr rQ +mgr sin 6 



dt 



t2 



d(m r 6r] + di mr'^O'^ Q 59 



. 



Suponiendo que 5r y 59 son ambas iguales a cero en ti y t2 la segunda integral 
es evidenteniente nula. Como 5r y 59 son completamente independientes, la 
primera integral puede ser cero solamente si 



.2 



m r —mr 9 —mg cos 9 + k{r — Tq) = 



mr 9 +2mr r0 +mgrsm9 = , 

pero estas son las ecuaciones de movimiento del sistema. 
b) Ejemplo de calculo de minimo 

Se trata de demostrar que la linea mas corta entre dos puntos cualesquiera 
pi J C sobre un cilindro es una helice. 

La longitud S de una li'nea cualquiera sobre un cilindro entre pi y p2 esta 
dada por 



P2 
PI 



1 + r^ 



d9 



V 



.1/2 



dz J 



dz , 



donde r, ^ y z son las coordenadas cili'ndricas usuales con r = cte. Puede 
determinarse una relacion entre 9 y z que le de a esta integral un valor 
extremo, por medio de 

d fd(t>\ d(t> 
Tz \W) ~ '09 

donde </>=[! + r^9''^] ^^"^ y 9' = f^, pero como d(l)/d9 = 0, 







f^(i + ^v^ 



"\v 



ci = cte. 



por esto r9' = C2. Y asi, r9 = C2Z + c^ que es la ecuacion de una helice. 
Supongamos que en pi se tiene 9 = y z = 0, entonces C3 = 0. En p2, 
liagase 9 = 92 y z = Z2 por tanto C2 = r92/z2, y r9 = {r92/z2)z es la 
ecuacion final. 
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2. MOVIMIENTO BAJO FUERZAS CENTRALES 

Prologo: Por rtiotivos astronomicos este fue el movimiento mas estudiado 
del punto de vista experimental en los primeros dos siglos de fisica moderna 
y representa un ejemplo principal para cualquier tipo de formalismo matem- 
atico. En su variante relativista, el problema de Kepler sigue siendo de gran 
interes. 



CONTENIDO: 

2.1 Problema de dos cuerpos: reduccion al problema de un cuerpo. 

2.2 Ecuaciones de movimiento. 

2.3 Ecuacion diferencial de la orbita. 

2.4 Problema de Kepler. 

2.5 Dispersion por un centro de fuerzas (con ejemplo). 
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2.1 Problema de 2 cuerpos: Reduccion al problema de un cuerpo. 

Considerenios un sistema monogeno de dos puntos niateriales de masas rrii y 
1712 ) en el cual las unicas fuerzas son las debidas a un potencial de interaccion 
V. Supondremos que V es funcion de cualquier vector entre ttt-i y m2, 
r2 — ri, o de sus velocidades relativas r2 — ri, o de derivadas superiores de 
r2 — ri. Este sistema tiene 6 grados de libertad y por lo tanto 6 coordenadas 
generalizadas independientes. 

Consideremos que estas son las coordenadas del vector de posicion del centro 
de masa, R, mas las tres componentes del vector diferencia r = r2 — ri. La 
Lagrangiana del sistema tendra entonces la forma: 

L = T{tl,r)-V{r,r,r, ). (1) 

La energia cinetica T es la suma de la energia cinetica del movimiento del 
centro de masa mas la energia cinetica del movimiento en torno al centro de 
masa, T: 



siendo 



T = -{mi+m2)K^ + T, 



TT 1 -2' I 1 -2' 

T= -miT^ + -m2r2. 



Aqui ri' y r2' son los vectores de posicion de las dos particulas relativas al 
centro de masa y estan relacionadas con r a travez de 

n = ■ r, r2 = ■ r , 2) 

mi + 771-2 ^T^l + ITl-2 



entonces T'toma la forma 



1 mim2 .2 
1 = r 

2 nil + ?n,2 



y la Lagrangiana total dada por la ecuacion (1) es: 

L=-[rm+ m2)R2 + ^— ^r^ - V{v, r, r, ) , (3) 

de donde definimos la masa reducida como 

17111712 , 1 1 1 

/i = ■ o - = — + — . 

mi + ?7l2 /U m-i 7TT-2 

Entonces nuestra ecuacion (^ se puede escribir como 

L = -{mi + m2)R^ + -/ur^ - V{r, r, r, ). 
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De esta ecuacion veinos que las coordenadas de R son ciclicas por lo que el 
centre de niasa estara fijo o se movera uniforniemente. 

Ahora, ninguna de las ecuaciones de movimiento para r contendra terminos 
donde aparesca R o R, este termino de la ecuacion es exactamente lo que 
tendrianios si tuvieramos un centro de fuerzas fijo (en el centro de niasa) 
con una particula a una distancia r de el con masa /U (masa reducida). 
Asi pues, el movimiento de dos cuerpos en torno a su centro de masa debido 
a una fuerza central se puede reducir siempre a un problema equivalente de 
un cuerpo. 

2.2 Ecuaciones de movimiento. 

Ahora nos limitaremos a fuerzas centrales conservativas, en donde el po- 
tencial es funcion solo de r, V{r), por lo que la fuerza solo estara dirijida 
a lo largo de r. Ya vimos que para resolver el problema solo necesitamos 
considerar el problema de una particula de masa m que se mueva en torno 
a un centro de fuerzas fijo, el cual sera el origen del sistema de coordenadas. 
Como el potencial solo depende de r, el problema tiene simetria esferica, 
es decir, cualquier rotacion, en torno a cualquier eje fijo, puede no tener 
efecto sobre la solucion. Por tanto, una coordenada angular que represente 
rotacion alrededor de un eje fijo debe de ser ciclica, lo cual da una simplifi- 
cacion considerable al problema. Debido a la simetria esferica el vector de 
momento angular total 

L = r X p 

se conserva. Se deduce, por tanto, que r es perpendicular a la direccion fija 
de L. Ahora si L = el movimiento debe ser a lo largo de una recta que 
pase por el centro de fuerzas, ya que para L = r y r son paralelas, cosa 
que solo se cumple en el movimiento rectilineo, por tanto el movimiento bajo 
fuerza central es siempre un movimiento piano. 

Ahora bien, tomando el eje z en direccion de L, el movimiento tendra sim- 
pre lugar en un piano normal al eje z. La coordenada esferica (j) tendra 
entonces el valor constante it/2 y podemos prescindir de ella en el estudio 
que sigue. La conservacion del momento cinetico proporciona 3 constantes 
de movimiento independientes. De hecho, dos de ellas, que expresan la 
direccion constante del momento cinetico, se han utilizado para reducir el 
problema de 3 grados de libertad a dos. La tercera corresponde a la conser- 
vacion del modulo de L. 
En coordenadas polares la Lagrangiana es 

L = -m{r'^ + r'^e^) - V{r) . (4) 
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como vimos 6 es coordenada ciclica cuya cantidad de momento canonico es 

el momento cinetico 

dL ^. 

Po = — r = mr t> , 
89 

entonces una de las dos ecuaciones de movimiento sera 

Pe = J^i^r^^) = ' (5) 

lo que nos conduce a 

mr 9 = 1 = cte , (6) 

donde / es la magnitud constante del momento cinetico. De la ecuacion (5) 
se deduce tambien que 

Se introduce el termino 1/2 por la razon de que {r'^9)/2 es la velocidad 
areolar (area barrida por el vector de posicion por unidad de tiempo) . 
La conservacion del momento cinetico es equivalente a decir que la veloci- 
dad areolar es constante. Tenemos aqui la demostracion de la segunda ley 
de Kepler del movimiento planetario: el radio vector barre areas iguales en 
tiempos iguales. Sin embargo debemos recalcar que la constancia de la ve- 
locidad areolar es una propiedad de movimiento debido a una fuerza central 
y no esta limitada a una ley de fuerza inversamente proporcional al cuadrado 
de la distancia. 
La ecuacion de Lagrange restante, para la coordenada r es 

0. (8) 

Designando por /(r) el valor de la fuerza, podemos escribir la ecuacion en 
la forma 

mr — mr9 = f{r) . (9) 

Utilizando la ecuacion (6), esta ecuacion se puede reescribir como 

mr 5- = f(r). (10) 

mr^ ^ ' ^ ' 

Basandonos en el teorema de la conservacion de la energi'a 

E = T + V = -mi'i-'^ + r'^9'^) + V{r) . (11) 

23 



d , , 


- mr9^ + -—- 


Jt^"^'^ - 


dr 



E es una constante de movimiento. Esto lo podemos deducir de las ecuaciones 
de movimiento. La ecuacion (10) la podemos escribir en la forma 



mr 



d 
dr 



1 /2 

2 mr'^ 



(12) 



Ahora multipliquemos por r ambos lados de la ecuacion 



d 1 d 

mrr = — {—mr) = 

dt^2 ' dt 



V{r) + - 



1 f 



o bien 



Por lo tanto 



d_ 

lit 



—mr + V{r) -\ — 



1 f 



2 mr"^ 



2 mr"^ 



. 



1 



-mr + V{r) + - 



1 f 



cte 



(13) 



2 ' ^ ' ' 2 mr"^ 

y ya que {P/2mr'^) = {mr'^6/2), la ecuacion (13) se reduce a (11). 

Ahora resolvamos las ecuaciones de movimiento para r y 9. Despejando r 

de la ecuacion (13), tenemos 




o bien 



dt 



V 



dr 



P 



2mr- 



r), 



2/ 2 



■^(E-V 

m \ 



«2 



(14) 



(15) 



2mr^ 



Sea ro el valor de r al timepo t = 0. La integral de los 2 miembros de la 
ecuacion toma la forma 



t 



dr 



TO 



21 2_ 
m 



{E-V 



2mr'^ ' 



(16) 



Esta ecuacion nos da t en funcion de r y de las constantes de integracion E, 
I y vq. No obstante se puede invertir, al menos formalmente, para dar r en 
funcion de t y de las constantes. Una vez hallada r, se deduce inmediata- 
mente 6 a partir de la ecuacion (6), que se puede escribir 



d9 



Idt 



mr^ 



(17) 
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Si 00 es el valor inicial de 6, entonces (17) sera 

^ = U ^2m+^o. (18) 

Asi pues hemos ya obtenido las ecuaciones de movimiento para las variables 
ry 9. 

2.3 Ecuacion diferencial de la orbita. 

Al tratar detalles concretos de problemas de fuerzas centrales reales conviene 
efectuar un cambio en la orientacion de nuestro trataniiento. Hasta ahora, 
resolver el problema significa liallar r y 6 en funcion del tiempo siendo E, 
I, etc. constantes de integracion. Pero muy a menudo, lo que realmente 
buscamos es la ecuacion de la orbita, es decir, la dependencia entre r y 6, 
eliminando el parametro t. En el caso de problemas de fuerzas centrales, esta 
eliminacion es particularmente sencilla ya que t solo figura en las ecuaciones 
de movimiento en forma de variable respecto a la cual se deriva. En verdad, 
una ecuacion de movimiento, (P), no hace sino darnos la una relacion definida 
entre una variacion infinitesimal dt y la variacion d9 correspondiente 

Idt = mr^d9. (19) 

La relacion correspondiente entre sus derivadas respecto a t y 9 es 

d I d 



dt mr"^ d9 



(20) 



Estas relaciones se pueden usar para convertir la ecuacion (|10|) en una 
ecuacion diferente para la orbita. Tambien se pueden solucionar las ecua- 
ciones de movimiento formalmente y llegar a la ecuacion de la orbita. De 
momento continuemos con la primera posibilidad. 

A partir de la ecuacion (pOl) podemos escribir la segunda derivada con re- 
specto a t 

(f_ _ d__l_ / d I 

dt"^ d0 mr"^ \d9 mr'^ 

y la ecuacion de Lagrange para r, (p!o|), queda en la forma 

I d f I dr\ I 



r"^ d9 \mr'^ d0 J mr^ 

Pero 

l_dr _ d{l/r) 

^d9~ d9~ 
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m . (21) 



Haciendo el cambio de variable u = 1/r, tenemos 



m 



dff^ 



+ u 



-f 



(22) 



Como 



la ecuacion (|2 



d dr d Id 

du dO dr v? dr 

puede escribirse en la forma 

d'^u 



de^ 



+ u 



m d / 1 
P du \u 



(23) 



Cualquiera de las dos ecuaciones (22) o (|23| ) es la ecuacion diferencial de la 
orbita si se conoce la fuerza / o el potencial V. Inversamente si conocemos 
la ecuacion de la orbita podemos seguir los pasos inversos y obtener f o V. 
Para una ley de fuerza particular cualquiera, la ecuacion de la orbita debe 
obtenerse integrando la ecuacion ( P2|) en una u otra forma. Puesto que ya 
se ha realizado la mayor parte del trabajo al resolver la ecuacion ([To|), solo 
queda eliminar t de la solucion ([T5|) por medio de (19), 



de 



Idr 



mr^ 



2/ J. 
m 



E - V{r) 



_P_ 



(24) 



dr 



ro 7-^ 



2 2/2mE _ 2mU _ 1 



~p W 



Haciendo el cambio de variable u = 1/r, 



du 



«0 .2/^ _ 2mt/ _ ^2 



(25) 



(26) 



lo que es la solucion formal para la ecuacion de la orbita. 

2.4 Problema de Kepler: Fuerza inversamente proporcional al 
cuadrado de la distancia 

La ley inversamente proporcional al cuadrado de la distancia es la mas im- 
portante de todas las leyes de fuerzas centrales por lo que le daremos un 
tratamiento detallado. En este caso la fuerza y el potencial son: 



/ 



k 



y y 



k 



(27) 
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Para integrar la ecuacion de la orbita sustituyamos (23) en (22), 

Hacemos el cainbio de variable y = u — ^ , para que la ecuacion diferencial 
quede en la forma 

cuya solucion es 

y = B cos{d - 9- ,) 

siendo B y 9' las correspondientes constantes de integracion. La solucin en 

fuencion de r es 

1 mk 

donde 

mk 
Podemos obtener la ecuacion de la orbita a partir de la solucion formal (|26 
A pesar de que este procedimiento es mas largo que resolver la ecuacion 
( psi) , resulta ilustrativo hacerlo ya que la constante de integracion e se evalua 
automaticamente en funcion de E j I. 
Escribamos ( P6|) en la forma 

= ff- I , '^'^ (30) 

J 2 2rnE 2mU „,2 

V~i^ T^ ^ 

donde ahora se trata de una integral indefinida. La cantidad 9' que aparece 



-^[l + ecos{9-9)], (29) 



en (30) es una constante de integracion determinada por las condiciones 
iniciales y no tiene por que ser el angulo inicial ^o al tiempo t = 0. La 
solucion a este tipo de integrales es 



dx 1 

arccos 



donde 



/3 + 27X 



(31) 



q = 0^ - Aa-f. 
Para aplicar este tipo de soluciones a la ecuacion (|30|) debemos hacer 

2m,E 2m,k 

Oi — p ) P — j2 ' "^ ~ ~ ' 
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y el discriminante q sera por lo tanto 

2 



(— ) V 



+ 



2Ef 



Con estas sustitucion 



mk"^ 
queda en la forma 



arccos 



L U 1 

mk 



1 + 



2EP_ 
mk-' 



Despejando n = 1/r, la ecuacion de la orbita resulta ser 
1 mk 



/2 



. 2 /, 2EP 
1+ WH nrCOSl 



m. 



A;2 



(32) 



Comparando ( ^2[) con la ecuacion (28) observamos que el valor de e es: 




(33) 



La naturaleza de la orbita depende del valor de e segiin el esquema siguiente: 



e > 1, E>0 

e = 1, ^ = 

e < 1, ^<0 

e = £; = - 



'W 



hiperbola, 
parabola, 
elipse, 
circunferencia. 



2.5 Dispersion por un centro de fuerzas. 

Desde un punto de vista historico, el interes acerca de las fuerzas centrales 
surgio en los problemas astronomicos del movimiento planetario. Sin em- 
bargo, no hay razon alguna para que solo las consideremos en este tipo 
de problemas. Otra cuestion que podemos estudiar mediante la Mecanica 
Clasica es la dispersion de particulas por campos de fuerzas centrales. Desde 
luego, si el tamario de las particulas es del orden atomico, debemos esperar 
que los resultados especificos de un tratamiento clasico sean a menudo in- 
correctos desde un punto de vista fisico, ya en que tales regiones suelen ser 
importantes los efectos cuanticos. A pesar de todo hay predicciones clasicas 
que siguen siendo validas con buena aproximacion. Mas importante aun, los 
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procedimientos de descripcion de los fenomenos de dispersion son los mis- 
mos en la Mecanica clasica que en la cuantica; podemos aprender a hablar 
el lenguaje igualmente bien basandonos en la Mecanica clasica. 
En su formulacion para un cuerpo, el problema de la dispersion se ocupa de 
la desviacion de particulas por un centra de fuerzas. Consideremos un haz 
uniforme de particulas -da igual que sean electrones, protones o planetas- to- 
das de igual masa y energia que inciden sobre un centro de fuerzas. Podemos 
suponer que la fuerza disminuye tendiendo a cero a grandes distancias. El 
haz incidente se caracteriza especificando su intensidad I (tambien llamada 
densidad de flujo), la cual da el niimero de particulas que atraviesan en 
unidad de tiempo la unidad de superficie colocada normalmente al haz. Al 
acercarse una particula al centro de fuerzas sera atraida o repelida y su orbita 
se desviara de la trayectoria rectilmea incidente. Despues de haber pasado 
el centro de fuerzas, la fuerza que se ejerce sobre la particula ira dismin- 
uyendo de manera que la orbita tendera de nuevo a tener forma rectilinea. 
En general, la direccion final del movimiento no coincide con la incidente 
y diremos que la particula se ha desviado o dispersado. Por definicion la 
seccion eficaz, (j{Q,), de dispersion en una direccion dada es 

a{n)dn = ^, (34) 

donde dN es el mimero de particulas dispersadas por unidad de tiempo en 
un elemento de angulo solido d^ en la direccion fi. A menudo, a a{i^) se le 
llama tambien seccion eficaz diferencial de dispersion. En el caso de fuerzas 
centrales debe haber una simetria total en torno al eje del haz incidente, por 
lo que el elemento de angulo solido podra escribirse en la forma 

dn = 2'Ksmede, (35) 

donde Q es el angulo que forman las direcciones desviadas e incidentes, al 
cual se le da el nombre de angulo de dispersion. 

Para una particula dada cualquiera, las constantes de la orbita y por lo 
tanto la magnitud de la dispersion, estan determinadas por su energia y su 
momento cinetico. Conviene expresar el momento cinetico en funcion de 
la energia y de una cantidad s llamada pardmetro de impacto que es, por 
definicion, la distancia del centro de fuerzas a la recta soporte de la velocidad 
incidente. Si uq es la velocidad incidente de la particula, tendremos 



I = muQS = s ■ v2rnE. (36) 
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Una vez fijadas E y s, queda deter minado univocamente el angulo de dis- 
persion Q. De momento supondrenios que valores diferentes de s no pueden 
llevar un mismo angulo de dispersion. Por tanto, el niimero de particulas 
dispersadas por un angulo solido dfl comprendido entre y G + d@ debera 
ser igual al niimero de particulas incidentes cuyo parametro de impacto este 
comprendido entre los valores correspondientes s y s + ds: 

2ttIs \ds\ = 27rcr(G)/ sin e\de\. (37) 

En la ecuacion (^) se han introducido los valores absolutos por que los 
niimeros de particulas tiene que ser siempre positivos, mientras que s y O 
varian a menudo en sentidos opuestos. Si consideramos s funcion de la 
energia y del angulo de dispersion correspondiente, 

s = s{e,E), 

la dependencia entre la seccion eficaz diferencial y vendra dada por 

ds 



a{Q) 



sin© 



dQ 



(38) 



A partir de la ecuacion de la orbita (|2^) se puede obtener directamente una 
expresion formal del angulo de dispersion. Tambien ahora, para mayor sen- 
cillez, consideraremos el caso de una dispersion puramente repulsiva. Como 
la orbita debe ser simetrica respecto a la direccion del periapside, el angulo 
de dispersion vendra dado por 

G = vr - 2* , (39) 

donde ^ es el angulo que forma la direccion de la asintota incidente con 
la direccion del periapside. A su vez, ^ puede obtenerse de la ecuacion 



(25) haciendo tq = oo cuando Oq = it (direccion incidente), por consiguiente 
6 = TT — ^ cuando r = Vm, distancia de mayor acercamiento. Facilmente se 
llega a 

* = / , . (40) 

Expresando / en funcion del parametro de impacto s (ec. (p6[)), resulta 

e = .-2r , ''' (41) 

: \^ _ Yill\ _ .2 
^ E 



30 



s^ 



o bien 



e 



vr 



sdu 



1 



v{u) 



s'^v? 



(42) 



Las ecuaciones (^) y ( ^^ rara vez se utilizan, a no ser en el calculo numerico 
directo del angulo de dispersion. No obstante, cuando se disponga de una 
expresion analitica para las orbitas, se puede a menudo obtener una relacion 
entre Q y s casi por simple inspeccion. 

EJEMPLO : 

Este ejemplo es historicamente muy importante. Se trata de la dispersion 
repulsiva de particulas cargadas por causa de un campo coulombiano. El 
campo de fuerzas dispersor es el creado por una carga fija —Ze al ejercerse 
sobre particulas incidentes que tienen carga —Z'e; por tanto, la fuerza se 
puede escribir en la forma 

es decir, se trata de una fuerza repulsiva inversamente proporcional al cuadrado 
de la distancia. Escribamos la constante 



k 



-ZZ'& 



(43) 



La energia E es mayor que cero y la orbita sera una hiperbola de excentri- 
cidad dada por 



'1 + 



2^/2 



/ 2Es \2 



1 + 



KzZ'e^ 






(44) 



m{ZZ'e^Y' 

donde hemos tenido en cuenta la ecuacion (|36D. Si se toma igual a vr el angulo 
^de la ecuacion (29), el periapside correspondera a = y la ecuacion de 
la orbita queda en la forma 



raZZ'e 
P 



[e cos 9-1]. 



(45) 



La direccion ^ de la asmtota de incidencia queda entonces determinada por 
la condicion r — > oo: 



1 



o sea, segiin la ecuacion (39), 



3c 


cos ^ = - 
e 


- - 


G 1 
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Luego 


















cot 


2©. 

2 ' 


= 6^-1, 


y utilizando la 


ecuacion 


([44) 














cot 


e _ 

2 ~ 


2Es 
ZZ'e^ ■ 



La relacion funcional buscada entre el parametro de impacto y el angulo de 
dispersion sera pues, 

de manera que efectuando la transformacion que exige la ecuacion (p3), 
encontramos que g(&) viene dada por 



La ecuacion (47) da la faniosa seccion eficaz de dispersion de Rutherford, 
quien la dedujo para la dispersion de particulas a por los nucleos atomicos. 
La mecanica cuantica da, en el limite no relativista, una seccion eficaz coin- 
cidente con este resultado clasico. 

En fisica atomica tiene mucha importancia el concepto de seccion eficaz total 
de dispersion ax cuya definicon es 

or = / a{n)d9. = 2tt [ a{e)de . 

JilT Jo 

No obstante, si intentamos calcular la seccion eficaz total para dispersion 
coulombiana sustituyendo la ecuacion (|4^) en esta definicion obtenemos un 
resultado infinito. La razon fisica de esto es facil de ver. Segiin su definicion, 
la seccion eficaz total es el niimero de particulas que, por unidad de inten- 
sidad incidente, se dispersan en todas direcciones. Aliora bien, el campo 
coulombiano constituye un ejemplo de fuerza de <Clargo alcance>>; sus efec- 
tos se extienden hasta el infinito. Las desviaciones muy pequeiias solo tienen 
lugar en el caso de particulas de parametro de impacto muy grande. Por 
tanto, todas las particulas de un haz incidente de extension lateral infinita 
se desviarian mas o menos y deben de incluirse en la seccion eficaz total de 
dispersion. Queda claro, pues, que el valor infinito de ax no es peculiar del 
campo coulombiano, tiene lugar en Mecanica clasica siempre que el campo 
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dispersor sea diferente de cero a todas las distancias independientemente de 
lo grande que sean. 
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3. CUERPO RIGIDO 



Prologo: Por las particularidades de su movimiento, el estudio del cuerpo 
rigido ha generado nuevas tecnicas y procedimientos matematicos intere- 
santes. 



CONTENIDO: 

3.1 Definicion. 

3.2 Grados de libertad. 

3.3 Tensor de inercia (con ejemplo). 

3.4 Momento angular. 

3.5 Ejes principales de inercia (con ejemplo). 

3.6 El teorema de los ejes paraleles (con 2 ejemplos). 

3.7 Dinamica del cuerpo rigido (con ejemplo). 

3.8 Trompo simetrico libre de torcas. 

3.9 Angulos de Euler. 

3.10 Trompo simetrico con un punto fijo. 
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3.1 Definicion. 

Un cuerpo rigido se define como un sistema de particulas cuyas distancias 
relativas estan obligadas a permanecer absolutamente fijas. 

3.2 Grades de libertad. 

Para describir el movimiento general de un solido rigido en el espacio tridi- 
mensional solo requerimos de 6 cantidades, por ejemplo: las 3 coordenadas 
del centro de masa medidas desde un sistema inercial y 3 angulos para especi- 
ficar la orientacion del solido (o de un sistema fijo en el solido con origen en 
el centro de masa) decimos que un cuerpo rigido en el espacio tiene 6 grados 
de libertad. 

El niimero de grados de libertad puede ser menor en los casos en que el 
solido esta sujeto a restricciones, por ejemplo: 

• Si el solido solo gira alrededor de un eje movil es de un grado de 
libertad (basta con un angulo). 

• Si el solido se mueve en el piano, su movimiento es mas general requiere 
de 5 cantidades (2 grados de libertad traslacional y 3 grados de libertad 
rotacional). 

3.3 Tensor de inercia. 

Consideremos que el cuerpo constituido por N particulas de masas m^, 
a = 1,2, 3..., N. Si el cuerpo rota con velocidad angular to alrededor de un 
punto fijo del cuerpo, y este punto a su vez se mueve a velocidad v respecto 
al sistema fijo (inercial), entonces la velocidad de la a -esima particula re- 
specto al sistema inercial esta dada por 

Vq, = v + w X ro,. (1) 

La energia cinetica de la a -esima particula es 

Ta = -rria^l , (2) 

Va = Va • Va = (v + (J X r^) • (v + 6J X r^) (3) 

= V • V + 2v • (w X r^,) + (uj X Ta) ■ {uJ X Ta) 

= ^r^ + 2v{L0XrJ + {uJXraf. (4) 
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donde 



Entonces la energia total es 

T = iMv2 + V [kjx J2a rriara] + I Ea "i„ (w X r«)2 . 
Si el origen esta fijo al solido lo elegimos en el centre de masas, entonces 



M 



por lo que 



1 



r=-Mv2 + -^m„(u;xrc 

a 
-'- — -^ trans ~r J rot 



donde 



1 



- trans 



^TTIqV^ 



1 



-Mv^ 



T, 



rot 



nll"^a 



CJ X To 



Ahora usaremos en la ecuacion (|8D la identidad vectorial 

(A X Bf = A^B^ - (A • B)2 
entonces la ecuacion nos queda en la forma 



r. 



ot 



^rua 



2 2 

w r 



(w • r„)^ 



(5) 
(6) 

(7) 
(8) 

(9) 



que en terminos de las componentes de a; y r 



1 



'-rot 



-.j: 



rrir 



k 



E^iXai ] J2^jX. 



ju^aj 



Ahora introducimos 



1 






T^rot = l^^^^ai UJiUJj5ij E xlk 
"^ a ij >- \ fc 



UJiUJjXaiXaj 



(10) 
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Trot = 2 ^ ^^^^ ^ "^" 
ij a 

Podemos escribir Tmt como 



Oij 2_i ^ak XoiiXdj 



fill 



1 



Tj-ot — -^ 2^ lij^i^j 



(12) 



donde 



E 



nir 



Oij 2_, ^ak ■^ai^aj 
k 



(13) 



Las 9 cantidades de lij constituyen las componentes de de una nueva can- 
tidad matematica que denotamos por {lij} y se llama Tensor de Inercia, 
{lij} se puede escribir convenientemente mediante un arreglo niatricial de 
(3 X 3) 

-^11 h2 hs 
{lij} = I hi I22 hz 

hi h2 hs 



Z2a "^a(^o2 + ^asJ 






2_,a ^^OiXalXaS 
Y,a''^a\^al + ^02) 



(14) 

Podemos notar que lij = Iji por lo tanto {lij} es un tensor simetrico, por 
lo tanto solo hay 6 terminos independientes. Los elementos diagonales de 
{lij} se llaman momentos de inercia con respecto a los ejes de coordenadas, 
los negativos de los elementos no diagonales se llaman productos de inercia. 
Para una distribucion continua de masa de densidad p(r), {lij} se escribe 
en lugar de (|l3|) como 



'«j 



P(r) 



Oij 2_i xj, XiXj 



dV. 



(15) 



EJEMPLO : 

Calcular los elementos hj del tensor de inercia {hj} para un cubo uniforme 

de lado 6, masa M, una esquina esta en el origen. 



b b b 



hi= I p 

V 



2 , 2,2 
Xi + X2 + X3 — XlXl 



dxidx2dx3 = p 



{x2 + X^)dxidX2dX3 . 



000 
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El resultado de la integral tres dimencional es In = ^{pb^)'^ = ^Mb'^. 

b b b 

h2 = / pi-xiX2)dV = -p / / ixiX2)dxidx2dx3 = -jpb^ = -jMb"^ 

V 

Vemos que las deinas integrales son las mismas por lo que 



L^ =ho = /■« = ^M62 



3 

J- 4 



entonces la matriz queda de la forma 



\Mb^ -\Mb^ -\Mh^ 



{li,} = I -\Mb' ^Mb' -\Mh' 
-\Mb'^ -iMb'^ iMb"^ 

3.4 Momento angular. 

El momento angular para el solido rigido constituido por A^ particulas rua 
esta dado por 

L =^r„ X Pa , (16) 

a 

donde 

Pa = maVoL = ma{uj X r^) . (17) 

Sustituyendo (|l7| ) en la ecuacion ([iq), tenemos que 

L =^mara X (w X Ta) ■ 
a 

Utilizando la identidad vectorial 

A X (B X A) = (A • A)B - (A • B)A = A^B - (A • B)A , 
tenemos 

a 

Tomando la z-esima componente del vector L 
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introduciendo la ecuacion 



J 



obtenemos 

^i = Y.a "^a [Jlj Sij^j Y.k ^ak) ~ [J^j Xaj^aj^j 
— l^a "^a l^j ^j"ij [Z-^k ^ak~-'^OiiXaj) 
= 2^j ^j 2^a "^a [Oij l^k ^ak~-^<^i'^cij) • 

Comparando con la ecuacion (|T3|) 

Li = 2_^ lijUJj . 



(18) 

(19) 
(20) 

(21) 



Esta ecuacion tambein se puede escribir en la forma 

L ={Iij}uj , 
o 



(22) 
(23) 



^ Li \ / 111 /i2 /l3 \ / ^1 

^2 = hi I22 hi, ^2 

\ L3 J y I31 /32 /as y y CJ3 

La energia cinetica rotacional, T^oi, se puede relacionar con el momento 
angular de la siguiente manera: multipliquemos la ecuacion ( ^l]) por ^uJi 



CJ. 



-Li — -LOi^L 



iUJ 



IJWJ , 



(24) 



sumando sobre todas las i da 



1 



2_^ ■::Li^i — 9 2^ hj^iOJj ■ 



2 2 



Al comparar esta ecuacion con (^), vemos que el segundo termino no es 
mas que T^ot, por lo tanto 



'^rot 



Y. n^i'^i = 9^ • 



UJ 



(25) 



Ahora, sustituimos (22) en la ecuacion (Eq), obtenemos una relacion entre 
la Trot y el tensor de inercia 



Trot = -UJ- {lij} -UJ- 



(26) 
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3.5 Ejes principales de inercia. 

Consideremos que el tensor de inercia {lij} es diagonal, es decir lij = lidij, 
la energia cinetica rotacional y el momento angular quedarian expresadas 
en la siguiente forma 

Trot = -^ 2^ hj^i^j 
1 



2 



7 ^ OijIiUJiOJj 



u 



^™t = jE^^^' (27) 



2 



y el momento angular 



Li — 2_^ lijUJj 
j 

= 2_^ ^ijIiiOj = liLOi 
J 

L = Icj. (28) 

Encontrar una expresion diagonal para {lij} equivale a encontrar un nuevo 
sistema de 3 ejes, en los cuales la energia cinetica y el momento angular se 
reducen a las expresiones ( p7| ) y (|2^), tales ejes se les llama Ejes Princi- 
pales de Inercia, es decir dado un cierto sistema inicial de coordenadas en 
el cuerpo, podemos pasar de el a los ejes principales mediante una transfor- 
macion ortogonal particular que, en consecuencia, se llama transformacion 
a los ejes principales. 



Igualando las componentes de (22) y (pS]), tenemos 



Li = lioi = hiLOi + I12UJ2 + /isi^s (29) 

L2 = IuJ2 = -^21^1 + -^^221^2 + -^231^3 (30) 

-^3 = I^'i = hi^i + h2^2 + hs^s , (31) 

las cuales son un conjunto de ecuaciones que se pueden reescribir 

(111 - I)uJi + I12UJ2 + /13W3 = (32) 

I21UJ1 + {I22 - I)UJ2 + -^23^3 = 

I^iuji + 132^2 + ihs - 1)^3 = . 
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Para obtener la solucion, el determinante del sistema debe ser cero 



(111 - I)uji I12LO2 -/'13W3 

l2i^i {I22 — 1)^2 hs^-i 

hl^l -^32^2 (-^33 — 1)^3 



. 



(33) 



El desarrollo de este determinante es un polinomio de grado 3 en I, 11a- 
mado polinomio caracteristico y la ecuacion ( |33D se llama ecuacion secular o 
ecuacion caracteristica. En la practica, los mementos principales de inercia, 
por ser los valores propios de I, se hallan buscando las rai'ces de la ecuacion 
secular. 

EJEMPLO : 

Determinar los ejes principales de inercia para el cubo del ejemplo anterior. 



Al sustituir los valores obtenidos en el ejemplo anterior en la ecuacion (|3 
obtenemos: 

({1/3-1) -i/3 -i/3 

-y (§/?-/) -1(3 I =0, 

V -i/3 -i(3 (1(3-1) 

donde /? = M6^, de donde obtenemos la caracteristica, 

entonces los eigenvalores o momentos principales de inercia son: 



h 



(3, h = h 



11 
12 



f3 



cuyos correspondientes eigenvalores son: 



1 

V 



Entonces la matriz que diagonaliza a {/jj} es: 

/i -f2 -fA 



1 



/ 
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{lij} diagonalizado es: 



U^,} 



diag 



(Aru..}A 



/ ¥ 











12 P 





V 





T^/3 



3.6 El teorema de los ejes paralelos. 

Supongamos que el sistema xi,X2,3;3 tiene su origan en el centre de masas 
del cuerpo rigido. Un segundo sistema Xi,X2,X3, tiene su origen en otra 
posicion diferente al sistema anterior, la unica condicion es que sean par- 
alelos, definamos los vectores r = (xi,X2,X3), R = {Xi,X2,X3) y a = 
(oi, a2, as), de tal manera que R = r + a o en terminos de sus componentes 



y\.i — Xi -\- di 



(34) 



Sean Jij las componentes del tensor de inercia respecto al sistema X1X2X3, 



'^ij ~ z^ 



rrir. 



k 



Sustituimos (34) en (^ 



'^ij ~ Z^ 



nir 



Oij /^[Xak ~r flfcj y^ai ~r CLijyXaj + CLj 

k 



(35) 



= [Ea "^a {Sij J2ki^ak)'^ " XaiXaj)] + J2a "^a {Sij J2k 4 " O-ittj) (36) 
+ Efc 2afc(5ij {J2a T^aXak) " O-j (Ea f^aXaj) " ^i (Ea ^^aXai)] ■ 

Pero la coordenada del centre de masa se define como 

y como habiamos dicho antes, el origen esta en el centro de masa 

(xi,X2,X3) = (0,0,0) . 

Ahora si tambien comparamos primer termino de ( |36[ ) con la ecuacion ([I^, 
tendremos: 

Jij = Iij + M{a^5ij- aittj) (37) 
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y entonces los elementos del tensor de inercia lij para el sistema del centre 
de masa estaran dadas por: 

lij = Jij — M{6ija — aittj) . (38) 

Este es el Teorema de los Ejes Paralelos. 

EJEMPLO : 

Calcular I^j para el cubo anterior respecto a un sistema paralelo al primer 

ejemplo y con origen en el centro de masa. 

Ya sabemos del ejemplo anterior que: 



/3 -y -y 



u^J} = I -y 1/5 -i/3 

-1/9 -i/5 §/? 



Ahora el vector a = (2,3)2) ^^^ ~ i^^' entonces usando la ecuacion (|3 
y el hecho que /3 = Mb"^ tenemos, 

111 = Jii - Mia"^ - al) = iMb"^ (39) 

I22 = J22 - M(a2 - ai) = iM62 (40) 

/33 = J33 - M(a2 - ai) = iM62 (41) 

/12 = J12 - M{-aia2) = (42) 

/12 = /i3 = ^23 = , (43) 



por lo tanto 



iM62 







{/} = I \Mb^ 




EJEMPLO : 

Consideremos el caso en el que el vector a = (0, g, 2) Y ^^ = T' entonces 
nuestro nuevo tensor de inercia seria: 

111 = Jii - M{a'^ - al) = (|M62) _ M (f - o) = \Mb'^ (44) 
I22 = J22 - M{a'' - al) = (|M62) _ M (§ - f) = ^Mb^ (45) 
/33 =J33-M(a2-ai) = (|M62)-M(f -f) = ^M62 (46) 

112 = Ji2-M(-aia2) = (-iM62) -M(0) = -iM62 (47) 
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Ii3 = Ji3 - M(-aia3) = (-iM62) - M(0) = -^Mb^ (48) 

/23 = J23 - M{-a2a3) = [-1Mb'') - M{\Mb'') = , (49) 

entonces {hj} es igual a: 



ih) 



( \Mb'' -iMb'^ -1Mb'' 

-1Mb'' ^Mb'' 

\^ -1Mb'' ^Mb' 



3.7 Dinamica del cuerpo rigido. 

La razon de cainbio respecto al tiempo del moinento angular L esta dado 
por: 

"'' / inercial 

Para la descripcion desde el sistema fijo al solido debenios usar la identidad 
operadora 

( — J = l-r] +'^^ • (^^) 

\ / inercial \ / cuerpo 



Aplicando este operador a la ecuacion (|50D 
fdL\ fdL\ 



4tl ,, \,tl +'^^^- '^^> 

/ inercial \ / cuerpo 



Entonces, en lugar de la ecuacion ( |50D tendremos 



a +.xL = N. (53) 

N / cuerpo 



Ahora proyectamos la ecuacion ( p3[ ) sobre los ejes principales de inercia, 
supongamos que estos son (xi, X2, X3), Trot y L se simplifican con tal eleccion, 
por ejemplo: 

Li = liCOi (54) 



la componente z-esima de (53) es 



—p + CijkUJjLk = Ni . (55) 
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Ahora proyectando sobre los ejes principales de inercia y utilizando la ecuacion 
(|5^, la ecuacion (|55| ) toma la forma: 

du>i 
li— + eijk^jujkh = Ni (56) 

ya que los elementos principales de inercia son independiantes del tiempo. 
Entonces, asi obtenemos un sistema de ecuaciones 

htOl + W2W3(/2 - -^3) = ^^1 (57) 

I2UJ2 + W3W1 (-^3 - h) = N2 
I3W3 + uJiuJ2{h - h) = N3 . 

Estas son las Uamadas Ecuaciones de Euler. 

EJEMPLO : 

Rodamiento y deslizamiento de una bola de billar. Demostrar que despues 

de un golpe horizontal la bola de desplaza resbalando una distancia 

_ 12-ug 

para despues empezar a rodar sin resbalar al tiempo 

t = — 
7At5' 

Al cesar la fuerza impulsiva las condiciones iniciales son: 

xo = 0, ±0 = Uq 

(/) = 0, = 0. 

La fuerza de friccion es 

entonces la ecuacion de movimiento es 

X = -ngM. (58) 

^=I,4> = Ns (59) 
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La ecuacion para L es 



donde Is es 



P(r) 



2 2 

X]^ — X2 



dxidx2dx'i = —Ma 



N3 = Ffa = ^iMga 
Sustituyendo en la ecuacion (|59|), tenemos 



2^. 



Ahora integrando las ecuaciones (58) y (pQ) una sola vez 

i; = — /i^t + Ci 

5 
«</* = 2/^5* + C2 

y aplicandoles las condiciones iniciales, tales ecuaciones quedan 

x{t) = —figt + Mo 

ant) = -fj-gt. 

Para que haya rodamiento puro, sin friccion se necesita que 

x{t) = a<j){t). 



(60) 

(61) 

(62) 
en la forma 

(63) 
(64) 

(65) 



De la ecuacion (64) y (pq) evaluadas en ti 

5 

-figh = -fj,gti + no 



2no 



h 



Ahora integrando de nuevo la ecuacion (p3l) y aplicando las 
iniciales tenemos que 



x(t) 



-^^9^ + uot . 

Evaluando la ecuacion ( pTJ ) y (p^) en el tiempo ti 

12^2 



(66) 
condiciones 

(67) 



X 



A9ng 
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. 5 

X = -Uo ■ 

3.8 Trompo simetrico libre de torcas. 

Un trompo simetrico es cualquier solido de revolucion. Si los mementos de 
inercia son 

Ii = h = h trompo esferico 

/i = /2 7^ I3 trompo simetrico 

Ii ^ I2 ^ I3 trompo asime'trico. 

Tomemos el caso del trompo simetrico Ii = I2 ^ h, en este caso el eje X3 
es el eje de simetria. Las ecuaciones de Euler proyectadas sobre los ejes 
principales de inercia son: 

/iL:;i+a;2u;3(/2-/3) = iVi (68) 

I2UJ2 + ^3^iih - h) = N2 (69) 

hu^ + uJiuj2{h-h)=N^- (70) 

Como el sistema que estamos considerando esta libre de torcas 

iVi = iV2 = iVs = , (71) 

utilizando el hecho que /i = I2 en la ecuacion ([7l|), obtenemos 

/iWi+a;2^3(-^2-/3) = (72) 

I2LO2 + uj-i'-^iih - /i) = (73) 

hd)^ = 0. (74) 

La ecuacion ( |7^ implica que 

61J3 = cte. 



Las ecuaciones (^) y (73) las reescribimos como: 



6ji = —Q.L02 donde = wa ( — j (75) 

L02 = — Ocui . (76) 
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Multiplicando la ecuacion ( [76| ) por i y sumandola a la ecuacion (^), tenemos 

((Ji + iL02) = -fl{uJ2 - ioJi) 
{uJl+iL02) = i^{uJi+iuj2)- 

Sea r]{t) = uJiit) + iCj2{i) entonces 

ri{t) - inr]{t) = , 



cuya solucion es 
Esto implica que 

entonces 



r]{t) = Aexp{i^t) 



{uJi + ioJ2) = ^ cos(Ot) + i sin(r2t) 



uji = Acos{9.t) (77) 

W2 = ^sm(m). (78) 

La magnitud del vector uj es 



cj = I Iwl I = \fuJi + W2 + 1^3 = Y ^^ + '^i = cte , 

esto significa que la magnitud de uo no canibia en el tiempo. Este vector 
realiza un moviniiento de precesion y la frecuencia de precesion esta dada 
por 

'h-Ii 



^ = t^3 , J 

ademas notamos que Vi es constante. 

Si le llamanios A al angulo entre u; y X^ las ecuaciones ( [77| ) y (78) toman la 

forma 

wi = usin Acos(ilt) 

iij2 = ^ sin A sin(r2t) 

oji = uj cos A , 

donde A = a; sin A. 

Para un solido de revolucion achatado Ii = I2 = I12 Y H > Ii- Para el caso 

de la tierra 

o f I3 — h2\ ^3 

i'm = ^3 
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Las observaciones indican un valor promedio de 14 meses ~ 450 dias. (Esto 
se debe a que no es un solido estrictamente, su estructura interna es li'quida). 

3.9 Angulos de Euler. 

Como sabemos una rotacion se puede representar por una matriz de rotacion 
A a traves de la ecuacion 

X = Ax . (79) 

X representa el conjunto de ejes del sistema rotado con respecto al sistema 
cuyos ejes represent amos por x'. La rotacion A puede llevarse a cabo por 
una sucesion de rotaciones "parciales" A = AiA2...A„. Existen muchas posi- 
bilidaddes para elegir estas A's. Una de ellas es el conjunto de angulos (p,6 
y ip llamados angulos de Euler, los cuales se generan a traves de la siguiente 
serie de rotaciones: 



• 



Una rotacion alrededor del eje X3 en un angulo if (en sentido positive). 
La matriz asociada a esta rotacion es: 



A(jO 



' cos (p sin if 
— sin Lf cos if 
^ 



• 



1 











COS 9 


sin 





— sin 6 


cos 



Una rotacion un angulo 6 alrededor del eje X'l (sentido positivo). La 
matriz asociada es: 



Xe 



Una rotacion en un angulo (p alrededor del eje X's (sentido positivo) 
la matriz asociada a esta rotacion es: 

/ coscj) sin(p 

— sin (f) cos (f) 

1 

La transformacion completa del sistema de ejes {X'l, X2, ^3} al sistema con 



ejes {Xi,X2,X3} esta dada por (|79|), donde 



^tf^^e^f 



49 



Haciendo el producto de matrices 
All = cos 9? cos 
A21 = 



donde 



cos o sm (p sm (p 



- sm 99 cos (p — cos y sm <p cos if 

A31 = sin sin (^ 

A12 = cos if sin (j) + cos 6 cos sin 93 

A22 = — sin 99 sin </> + cos cos (j) sin 99 

A32 = — sin Lp cos (p 

A13 = sin 99 cos 

A21 = cos if sin 

A33 = cos 9 



A 



All A12 Ai3 
A21 A22 A23 
A31 A32 A33 



Ahora, consideremos el hecho de que: 

• <j) esta dirigida a lo largo del eje X'^ (fijo). 

• 6 esta dirigido a lo largo de la linea de nodo. 

• if esta dirigido a lo largo del eje X3 (cuerpo). 



Podemos escribir 3 componentes de cada uno de los 3 vectores en el sistema 
{Xi,X2,X{\ como: 

(pi = (p sin 6 sin 99, 9i = 6 cos 99 991 = 
(p2 = (/> sin cos 99, 02 = —9 simp 992 = 
<pi = (pcose, 03 = 03 = 0, 

entonces 

LO = (p + 6 + if 

h+Oi + pi) , (<p2 + O2 + P2) , (<p3 + O3 + ifs) . 
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Entonces las coniponentes de uj son: 

uji = (j) sin 9 sin ip + 9 cos 99 
UJ2 = (f) sin 9 cos c^ — sin c^ 
1^3 = (t)cos9 + ip. 

3.10 Trompo simetrico con un punto fijo. 

Como ejemplo mas complicado de la aplicacion de los metodos de la dinaniica 
del cuerpo rigido, vamos a considerar el movimiento de un cuerpo simetrico 
en un campo gravitacional uniforme cuando un punto del eje de simetria 
este fijo en el espacio. 

El eje de simetria es, desde luego, uno de los ejes principales y lo tomaremos 
como el eje z del sistema de coordenadas solidario al cuerpo. Como hay 
un punto fijo, la configuracion del trompo quedara determinada por los tres 
angulos de Euler: 9 da la inclinacion del eje z respecto a la vertical, (j) mide 
el acimut del trompo respecto a la vertical, mientras que ^p es el angulo de 
rotacion del trompo respecto a su propio eje z. La distancia del centro de 
gravedad al punto fijo sera representada por /. para obtener una solucion 
del movimiento del trompo vamos a utilizar el metodo de Lagrange en vez 
de las ecuaciones de Euler. 
La energia cinetica es: 

o bien, en funcion de los angulos de Euler 

T = i/i(02 sin^ 9 + 9^) + ^/3((/>cos 9 + ^f , 

donde ban desaparecido los terminos en los que figuraban uj\ y wl- Conoce- 
mos un teorema elemental segiin el cual en un campo gravitatorio constante 
la energia potencial es la misma que se tendria si el cuerpo estuviera con- 
centrado en su centro de masa, pero vamos a dar una demostracion formal 
del mismo. La energia potencial del cuerpo es la suma extendida de todas 
sus particulas: 

V = -rriiYi ■ g , (80) 

donde g es el vector constante que representa la aceleracion de la gravedad, 
segiin como se define el centro de masa, esto es equivalente a 

V = -MR, • g, (81) 
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lo que deinuestra el teorema. La energia potencial en funcion de los angulos 
de Euler es: 

V = Mgl cos e, (82) 

con lo que la Lagrangiana sera 

L = ]-Ii{^'^ sm^ e + e'^) + ]-h[^ cos e + ^f -Mgl cose. (83) 

Notamos que y C/3 son coordenadas ciclicas, por lo tanto p,^ y p^p son 
constantes del movimiento. 

dL 

Pip = -^ = hiif + <j) cos 9) = cte (84) 

y 

O T 

P(j) = — - = Ii(j)s\T? 6 + /3((/)cos^ 9 + If COS 9) = cte. (85) 



De la ecuacion (84) despejamos (p 



p^-hcj) cos 9 
LP = -^ , (86) 

sustituimos en la ecuacion ( ^5| ) 

^^ T 1 ■ 20 , T fl 2/) , Pv-h4'Cos9 
Pff) = — - = liqjsui 9 + l2,[(p cos a -\ — cos 9) = cte 

d(j) h 

Ptj, = Ii(j>sm'^ 9 + p^cos9 , 

de donde obtenemos 



Py, cos 9 



(87) 



h sin^ 9 
Sustituyendo en (|86|) 

P-fi P<h - Pu, COS 9 

V? = -7- ^"o COS 9. 

h h sin^ 9 

Ahora, como el sistema es conservativo, otra constante de movimiento es la 
energia. 

E = T + V = ]-Ii{^^ sin^ 9 + 9'^) + ]-h{4>cos9 + ^pf + Mgl cos 9. 
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La cantidad I3UJ3 = p^p es una constante de moviento, niultiplicando esta 
constante por Pipco^ se obtiene 



7-2, ,3 



p1^3 
p1^3 

2h 



La cantidad 2-^3'^! ^^ una constante, entonces podemos definir la cantidad 



E' 



de donde identificamos 



1 2 
E - 2-^3^3 = cte 

-hO'^ + -/i(/>2 sin^ 61 + Mc/Z cos 9 , 



V{e) = -0^ sin^ 6 + M5/ cos ( 



^ ^ 2 ' V /i sin2 



sin^6' + Msr/cos6'. 



Entonces i?'es: 



a2 



E-=-hd' + V{d) 



De esta ecuacion despejamos 9 e integranios para obtener 



(89) 



h 



de donde obtenemos 



m 



(E'-vie)) 



de 



de 
IE ' 



^{j^){E-v{e)) 



(90) 



Al realizar la integral de la ecuacion (^) se obtiene t = f{0), de donde en 
principio podemos despejar y obtener 9{t). Entonces 6{t) se sustituye en 
las ecuaciones para </> y 9? (ecs. (q7|) y (H)) y al integrarlas obtenemos la 
solucion completa a nuestro problema. 
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4. OSCILACIONES PEQUENAS. 



Prologo: Una forma muy cormin de movimienteo en los sistemas niecanicos, 
son las pequenas oscilaciones. Estas las encontramos en sistemas tales 
como vibraciones atomicas, moleculares, circuitos electricos, acustica. Todo 
movimiento alrededor de las posiciones de equilibrio estable, es el llamado 
vibratorio. 



CONTENIDO: 

4.1 OSCILADOR ARMONICO SIMPLE. 

4.2 OSCILADOR ARMONICO FORZADO. 

4.3 OSCILADORES ARMONICOS AMORTIGUADOS. 

4.4 MODOS NORMALES. 

4.5 RESONANCIA PARAMETRICA. 
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4.1 OSCILADOR ARMONICO SIMPLE. 

Un sistema se encuentra en equilibrio estable cuando su energia potencial 
U{q) es minima; al separarlo de esta posicion se origina una fuerza —dU/dq 
que tiende a devolver al sistema al equilibrio. Sea go el valor de la coorde- 
nada generalizada correspondiente a la posicion de equilibrio. Al desarrollar 
U{q) — U{qo) en serie de potencias de Taylor de q — qo para pequeiias desvia- 
ciones de la posicion de equilibrio 

U{q)-U{qo)^h{q-qof , 
donde: 

— - n 
dq 

U{q) = 0, 

es decir: no hay fuerzas externas que actuan sobre el sistema y se ha escogido 
el nivel de referenda de tal modo que coincide con la posicion de equilibrio; 
ademas de despreciar terminos de orden superior. El coeficiente k representa 
el valor de la segunda derivada U(q) para q=qo. Por simplificacion haremos 
la siguiente designacion 

x = q-qo 
con lo que la ecuacion de energia potencial toma la forma: 

U{x) = -kx^ . (1) 

La energia cinetica de un sistema es en general de la forma 

1 -2 

T = -mx (2) 

2 ^ ^ 

con (|l]) y (y) obtemos la expresion para la Lagrangiana de un sistema que re- 
aliza oscilaciones lineales (a tal sistema se le llama frecuentemente oscilador 
lineal) : 

1-212 

L = —m X kx . (3) 

2 2 ^ ' 

La ecuacion de movimiento correspondiente a esta L es: 
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m X +kx = , 
o bien 

i +w^x = , (4) 

donde w'^ = ^/kjm. La ecuacion diferencial tiene dos soluciones independi- 
entes: coswt y senwt, asi, formamos la solucion general: 

X = c\ cos wt + C2sent(;t , (5) 

o bien, podemos expresar la solucion de la forma: 

x = acos(wt + a) . (6) 

Puesto que cos(wi + a) = cosu'icosa — senwtsena, la comparacion con 
(P) muestra que las constantes arbitrarias a y a estan relacionadas con los 
coeficientes c\ y C2 de la forma: 



(cf + C2), y tana = — ci/c2 . 

Asi: un sistema en las proximidades de su posicion de equilibrio estable, 
ejecuta un movimiento oscilatorio armonico. El coeficiente a en (^) es la 
amplitud de las oscilaciones, y el argumento del coseno su fase; a es el 
valor inicial de la fase, y depende evidentemente de la eleccion del origen de 
tiempos. La magnitud w es la frecuencia angular de las oscilaciones, esta no 
depende de las condiciones inicales del sistema por lo cual es la caracteristica 
fundamental de las oscilaciones. 
A menudo la solucion es expresada de la forma 

X = re [^exp(^^i;t)] 

donde A es la amplitud compleja, su modulo es la amplitud ordinaria: 

A = oexp(za) . 

La energia de un sistema que realiza pequefias oscilaciones es: 

^ 1 -2 1, 2 

E = -m X +-kx , 
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o sustituyendo (S) 



2 2 

E = —mw a . 
2 



Ahora, considereinos el caso para un niimero n de grades de libertad. En 
este caso seguiremos considerando que la suma de las fuerzas que actuan es 
cero, con lo que 

dU 
Qi = -j- = 0. (7) 

dqi 

Procediendo de igual forma que en el caso de un solo grado de libertad, re- 
alizamos una expancion en series de Taylor para la energia potencial donde 
ahora consideramos un mmimo para qi = QiQ. Introduciendo pequenos de- 
splazamientos 

Xi = Qi- QiO , 
al realizar la expancion en series 



(8) 
Bajo las mismas consideraciones que se hicieron en (||), llegamos a la sigu- 
iente relacion: 

U{qi,q2, ..., qn) = U = -Y^ kijXiXj . (9) 

De la ecuacion @ se puede notar que kij = kji, es decir son simetricos 
con respecto a sus indices. Ahora considerenios la situacion para la energia 
cinetica. Esta en general es de la forma 

—ciij{q) XiXj , 

donde las aij solo son funciones de las coordenadas. Al designar estas Ojj = 
rriij la energia cinetica es de la forma 
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Una vez conociendo las energias tenemos que la Lagrangiana para sistemas 
con n grades de libertad es de la forma 






fill 



Esta Lagrangiana lleva a las ecuaciones diferenciales de movimientos siniul- 
taneas 



d dL dL 



o bien 



dt Q X 9xi 



/ J V^'i^j *^i ' ij'^j ) — ^ • 



(12) 



(13) 



Tenemos asi un sistema de ecuaciones diferenciales lineales y homogeneas, 
las cuales pueden pueden ser considerasdas como las n componenetes de la 
ecuacion matricial 



{M){X) + {K){X)=0 
donde las matrices estan definidas por: 



(M) 



/ mil "1-12 ••• min \ 

17121 W-22 ••• '"T'2n 



(K) 



nir 



\ nini m„2 ••• 

/ kii ki2 ... hn \ 
k21 k22 ... k2n 



\ kfil kn2 ... knn / 

f XI \ 
d^ 



J 



(X) 



dt^ 



X2 



\Xn ) 



(14) 



(15) 



(16) 



(17) 
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(X) 



X2 



(18) 



\ Xn J 

De manera analoga al sistema con un grado de libertad, buscamos n fun- 
ciones incognitas Xj{t) de la forma 

Xj =Ajexp{iwt) , (19) 

siendo A^ constantes a determinar. Sustituyendo (^) en (^) y dividiendo 
todo entre exp{iwt), se obtiene un sistema de ecuaciones algebraicas lineales 
y homogeneas, a las que deben satisfacer Aj. 



^{-w'^niik + kik)Ak = 



(20) 



Para que este sistema tenga soluciones distintas de cero, el determinante de 
sus coeficientes debe anularse. 



"u 



w rrii 



. 



(21) 



Esta es la ecuacion caracteristica y es de grado n con respecto a w'^. En 
general, tiene n raices distintas reales y positivas Wa (a = 1,2, ...,n). Las 
magnitudes Wa se llaman frecuencias propias del sistema. Multiplicando por 
A* y sumando sobre i se tiene 

^(--w^my + kij)A*Aj = , 



de donde 



w 



/ ^ kijA^ Ail y ^ vfiijAi Ai . 



Como los coeficientes kij y rriij son reales y simetricos, las formas cuadraticas 
del numerador y denominador de esta expresion son reales, y al ser esencial- 
mente positivas, 'uP' es igualmente positivo. 

EJEMPLO 
Como ejemplo se modelaran las ecuaciones del movimiento del pendulo 
doble. La energia potencial para este es (el problema posee dos grados 
de libertad) 
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U = migli(l — cos^i) + m2gh{l — cos^i) + ^,25/2(1 — 003^2) 
Al aplicar la expancion (P), se tiene 



U = -(mi + m2)gh9l + -m2gl29l . 



Al comparar con 



identificamos 




kii = 


(mi +m2)/i 


ki2 = 


^21=0 


k22 = 


m2gl2 ■ 



Para la energia cinetica se encontro 

1 -2 1 -2 

r = 2("^i + "^2)/? 6^1 +-m2ll 02 +m2lil2 OiOi ■ 

Identificando terminos al comparar con ([Tol ) 

mil = (mi+m2)/i 
fni2 = m2i = m2lil2 
m22 = m2ll. 

Al sustituir las ecuaciones de las energias en (|ll|) se obtiene la Lagrangiana 
para el oscilador de pendulo doble y como resultado final: 

/mil mi2\('ei\^(kii W ^1 \ 
1^ m2i m22 ) \'e2 J V k22 ) \02 ) 

Estas son las ecuaciones de movimiento en este caso. 



4.2 OSCILADOR ARMONICO FORZADO. 

Si un sistema oscilatorio se somete a la accion de un campo externo vari- 
able, son las llamadas oscilaciones forzadas. Como consideramos pequeiias 
oscilaciones, entonces esperamos que la accion del campo exterior sea debil. 
Ademas de su energia potencial propia, el sistema posee en este caso una en- 
ergia potencial Ue{x,t) debida al campo exterior. Desarrollando esta ultima 
en serie de potencias de la pequeiia magnitud x: 
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Ue{x,t)^UeiO,t)+X 



dUe 

dx 



x=0 



El segundo termino, es la fueza exterior que actiia sobre el sistema en su 
posicion de equilibrio, desisgnado este termino como F{t). Con esto, la 
Lagrangiana para este sistema es de la forma 



T 1 -2 1, 2 

L = —m X ' — 

2 



-kx'^ + xF{t) 



(22) 



La ecuacion de movimiento correspondiente es 

m X +kx = F{t) , 
o bien 

X +w'^x = F{t)/m 



(23) 



donde w es la frecuencia para las oscilaciones libres. La solucion general a 
esta ecuacion es de la forma 

X — X}^ -f- Xp , 

es decir, de una parte homogenea y una solucion correspondiente a un caso 
particular. Analizando el caso para el cual la fuerza exterior es funcion 
periodica simple del tiempo, de frecuencia 7 de la forma 

F(t) = /cos(7t + /3) . 
Al hallar la integral particular para la ecuacion |2^ en la forma xi =h cos(7t+ 



j3) y al sustituir, se tiene que b = f/m{w'^ 
soluciones, tenemos que la solucion total es 



7 ), que al juntar ambas 



X = a cos{wt + a) + f/m{w —7) cos(7t + /3) . 



(24) 



El resultado muestra una suma de dos oscilaciones: una debida a la frecuen- 
cia propia y otra con la frecuencia de la fuerza exterior. 
La ecuacion (^) puede ser integrada en forma general para una fuerza ex- 
terior arbitraria. Escribiendo la ecuacion de la forma 

— ix +iwx) — iwix +iwx) = — F(t) , 
at m 
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haciendo ^ =x +iwx, se tiene 

— ^ — iw$, = F{t)/m . 

La solucion a esta ultima es del tipo ^ = A(t) exp{iwt); para la funcion A{t) 
se obtiene 

A= F{t)exp(—i'wt)/m . 

Al integrarla se obtiene la solucion 

/■* 1 

C = exp(iwt) —F(t)exjp(-iwt)dt + Co- (25) 

Jo rn 

Esta es la solucion general buscada; la funcion x(t) esta dada por la parte 
imaginaria de esta ultima dividiendo por w. 

EJEMPLO 
Como empleo de la ecuacion anterior se muestra el siguiente ejemplo. 
Determinar la amplitud final de las oscilaciones de un sistema bajo la accion 
de una fuerza exterior tal que Fq = cte durante un tiempo limitado T. Para 
este intervalo de tiempo se tiene 



^ = — exp(iti;t) / eyip{—iwt)dt 
m Jo 



Fo 

[1 — exp(— itut)] exp{iwt) 



iwm 
y el modulo cuadrado da la amplitud. De la relacion |^| = a^iv^ con lo cual 

2Fn 



a 



'^0 A rr\ 

-^sen[-wl ) 



mw"^ 2 

4.3 OSCILADOR ARMONICO AMORTIGUADO. 

En las secciones anteriores consideramos solo la carencia o la presencia de 
fuerzas externas, para los casos del oscilador armonico simple y forzado; 
respectivamente. Es decir, el movimiento tenia lugar en el vacio o bien, que 
la influencia del medio en el movimiento era despreciable. En la realidad, 
cuando un sistema se mueve a traves de un medio, este ofrece resistencia que 
tiende a retardar el movimiento. La energia del sistema se disipa (ya sea en 
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forma de calor 6 de algun otra forma de energia). Primero analizamos como 
afecta este fenomeno a las oscilaciones simples. 

El modo en que este medio afecta al movimiento es por fuerzas de roza- 
miento. Si esta fuerza disipativa es lo suficientemente pequeiia, podemos 
desarrollarla en potencias de la velocidad. El termino de orden cero del de- 
sarrollo es nulo, ya que ninguna fuerza de rozamiento actiia sobre un cuerpo 
enreposo, por lo que el primer termino que no se anula es proporcional a la 
velocidad, ademas despreciando terminos de orden superior. 

fr = -ax , 

donde x es la coordenada generalizada y a un coeficiente positivo; el signo 
menos indica que es en sentido opuesto al movimiento. Aiiadiendo esta 
fuerza a la ecuacion de movimiento 

m x= —kx — ax , 

o bien 

x= —kx/m — ax /m . (26) 

Haciendo k/m = w"^ y a/m = 2A; donde Wq es la frecuencia de las oscila- 
ciones libres del sistema y A es el coeficiente de amortiguamiento. Con lo 
anterior 

X +2A X +wlx = . 

La solucion para la ecuacion anterior es de la forma x = exp(ri); al sustituir 
esta en la ecuacion anterior obtenemos la ecuacion caracteristica para r. De 
tal modo 

r^ + 2A + u-^ = , 

de donde 



ri,2 = -A±^(A2-u;2), 
con lo que la solucion general a la ecuacion de movimiento es 

X = ci exp(rit) + C2 exp(r2t) . 
De las raices de r podemos considerar los siguientes casos especiales: 
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(i) X < Wo- Se tienen raices imaginarias conjugadas. Con lo que la 
solucion es 



X = re < Aexp 



-At + i^K - A2) 



siendo A una constante compleja arbitraria. La solucion puede ser escrita 
en la forma 



a exp(— At) cos{wt + a) siendo 



w 



VK - A2) , 



(27) 



donde a y a son constantes reales. De este modo, puede decirse que una 
oscilacion amortiguada es como una oscilacion armonica cuya amplitud de- 
crece exponenecialmente. La rapidez de disminucion de la amplitud esta 
deter minada por el exponente A y la frecuencia w es menor que las oscila- 
ciones libres en ausencia de rozamiento. 

(ii) X > Wo- Entonces los dos valores de r son reales y negativos. La 
forma general de la solucion es: 



X = ci exp 



X-^{X^-w^o) 



t> + C2 exp 



A + v'(A2-^2) 



Si el rozamiento es muy grande, el movimiento consiste en una disminucion 
monotona, que tiende asintoticamente (cuando t — > cxd) a la posicion de 
equilibrio (sin oscilacion). Este tipo de movimiento se llama aperiodico. 
(iii) A = Wo- Se tiene que r = —A, cuya solucion general es de la forma 

X = (ci + C2t) exp(— At) . 

Si generalizamos para sistemas con n grados de libertad, las fuerzas de roza- 
miento generalizadas correspondientes a las coordenadas Xi son funciones 
lineales de las velocidades 



Jr,i — Z^ 



Oiij X{ , 



(28) 



con aik = aki se puede escribir como 



fr. 



dF 

dxi 



donde F = ^^ J2i,j o^ij ^i^j Y se le llama funcion disipativa. La ecuacion 
diferencial se obtiene al sumar estas fuerzas a la ecuacion (O ) 
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/ A'^ii ^j '^'^ij^jj — / . 



CXij Xi . 



(29) 



Haciendo en estas ecuaciones 



Xk = Ak exp(rt) 



y al sustituir esta ultima en (|29| ) y dividiendo por exp(rt), se tiene el siguiente 
sistema de ecuaciones algebraicas lineales para las constantes Aj 

^{rriijr'^ + Oj^r + %)Aj = . 



Igualando a cero el determinante de este sistema, se encuentra la ecuacion 
caracteristica para este sistema. 



rriijr + aijr + kij 
Esta es una ecuacion en r de grado 2n. 







(30) 



4.4 MODOS NORMALES. 

Antes de definir los modos normales, reescribiremos la ecuacion (^) de la 
siguiente manera 



X) + K\X) = 



M 



donde \X) es el vector n-dimensional cuya matrix de representacion es ([iq); 
M y K son dos operadores que tienen la representacion matricial definidas 
por (15) y (16) respectivamente. La ecuacion antes descrita es una ecuacion 
con operadores. Dado que M es un operador no sigular y simetrico para el 
cual el operador inverso M~^ y el operador M^'^ y M"^'^ existen. Con lo 
anterior podemos expresar la ecuacion de operadores en la forma 






M^/2 IX) = -M-^^^KM-^/^M^^^ \X) 



o de la forma mas compacta 



dt2 



X 



X 



(31) 
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donde 



Como M^^' ^ y K son operadores simetricos, entonces A es igualmente simetrico. 
Si empleamos eigenvectores ortogonales como base vectorial (por ejemplo el 
espacio tridimencional), la representacion matricial del operador puede ser 
diagonal de la forma 

Trataremos el siguiente problemas de eigenvalores 

\\p^) = Xi\pi) , (32) 

donde \pi) representa un conjunto de eigenvectores mutuamente ortogonales; 
o bien 

M-^/^KM'^'^ \pi) = A, \pi) . 

Los eigenvalores son obtenidos al multiplicar ambos lados por {pi\, con lo 
cual 

^ _ {pi\M-^/^KM-y^\p,) 
M — 



{pi \pi) 

Dado que las energias potencial y cinetica son cantidades positivas, se tine 
que 

{p,\M~^/^KM'^/^\pi))0 
y por lo tanto 

Xi>0. 
Esto nos permite el conjunto 

Aj = wf . 
Si expresamos el vector x) en terminos de estos eigenvectores de A, 
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-^ 



J2'J'\^ 



donde 



y^ = {pi \x) . (33) 

Al insertar este resultado en la ecuacion de movimiento (pll), se tiene 






-A 



X 



Ea 



iVi \Pi) 



El producto escalar de esta ecuacion con el eigenvector constante {pj\, pro- 
duce la ecuacion de movimiento para las coordenadas generalizadas yj 






-WjVj 



La solucion para esta ecuacion es de la forma 

Hj = Aj cos{wjt + (j)j) . (34) 

En base a estas nuevas coordenadas para el movimiento armonico de un sis- 
tema de particulas, pueden obtenerse un conjunto de ecuaciones de movimiento 
generalizadas e independientes. La relacion entre estas y/s y Xi 's esta dada 
por (ID 

Vj = Pjl ^1 +Pj2 ^2 +••• + Pjn Xn ■ 

Los componentes pji (/ = 1, 2, .., n) son determinados al resolver el problema 
de eigenvalores de la ecuacion (p2|). Las nuevas coordenadas son referidas 
como las coordenadas normales, y las WjS como las frecuencias normales. 
La forma equivalente (^4|) en matrices es 



/ xf \ 

U) 

2 



XI 



\ 



?i3) 



Aj cos{'Wjt + (pj) 



( Pn \ 

Pj2 
V Pjn J 



(35) 



Esta ultima refiere los modos normales de vibracion del sistema. Una de las 
razones de haber introducido coordenadas yj's se aprecia al observar que la 



expresion para la energia cinetica no cambia si se rotaran los ejes del nuevo 
sistema no cambian 

1 



^=oE^.-y. • 



i=i 



EJEMPLO 



Suponiendo un arreglo niatricial de la forma mostrada, la cual representa 
las ecuaciones de movimiento obtenidas para un sistema 




5 1 \ 


(x. 


2 


X2 


10 5/ 


\ X3 



Comparando con (pl|), identificamos el valor del operador A. Al encontrar 
los egienvectores empleamos (|3^ con lo cual 



5 1 \ 


( ^' 


2 


P2 


10 5, 


V P3 




Determinamos la ecuacion caracteristica para las Aj. De manera que 

det(A - XJ) = , 
al sustituir valores 



5-A 1 

2-A 

1 5-A 



. 



Al resolver la ecuacion se tiene que Aj = 2, 4, 6. Para A = 4 



/5 l\ 


( PI 


2 


P2 


[l 5] 


\P3 




se tiene el siguiente conjunto de ecuaciones 
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(5-4)pi+p3 = 

2/92 - 4/92 = 

/9i + (5 - 4)p3 = . 
Teniendo en cuenta la condicion de normalizacion, se obtienen los valores 



Con lo anterior 



Pi 

P2 

|PA=4) 



De igual modo se obtienen 



-P3 



1 

V2 



V2 



\PX=6) 



\PX=2) 




Asi, el nuevo espacio vectorial esta determinado por 

I V2 V2 

\pi 



1 

V2 


1 

V2 



1 

V2 



1 

V2 



de donde 




(Pi 



Asi, las coordenadas normales son dadas por (33) 

Vol 
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4.5 RESONANCIA PARAMETRICA. 

El importante fenomeno de resonancia parametrica se presenta al tener un 
sistema que se se encuentra en un estado de reposo (en la posicion de equi- 
librio X = 0) y es inestable; es decir, bastara una separacion de esta posicion 
por pequeiia que sea para provocar un desplazamiento x rapidamente cre- 
ciente con el tiempo. Se diferencia de las resonancias ordinarias, en las cuales 
el desplazamiento crece con el tiempo (proporcional a t). 
Los parametros de un sistema lineal son los coeficientes TTi y A: de la La- 
grangiana (^); si estos son funcion del tiempo, la ecuacion del movimiento 
es: 

— {mx) + kx = Q. (36) 

Si consideramos la masa constante, la ecuacion anterior toma la forma 

-7-^+w^(t)x = . (37) 

La forma de la funcion w{t) esta dada por las condiciones del problema; 
suponiendo que la funcion es periodica de frecuencia 7 (y de periodo T = 
27r/7). Lo que significa 

w{t + T) = w{t) , 



por lo cual, toda ecuacion del tipo ( |37D es invariante con respecto a la trans- 
formacion t — > t+T. Asi, si x{t) es una solucion de esta, la funcion x{t+T) es 
tambien solucion. Con lo anterior, sean xi{t) y X2{t) dos integrales indepen- 
dientes de (^), estas deben transformarse en ellas mismas en combinacion 
lineal cuando se sustituye t ^ t + T. La forma de ello es 



xi{t + T) = fiix{t) (38) 



o en forma general 



x,{t) = f,fF{t) 
X2{t) = fA^^G{t) 
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donde F{t) y G{t) son funciones puramente periodicas del tiempo (de peri'odo 
T) . La relacion entre estas constantes se obtiene al nianipular las ecuaciones 
siguientes 

xi +w'^{t)xi = 

X2 +W^{t)x2 = . 

Multiplicando por X2 y xi respectivamente, y restandolas miembro a miem- 
bro, se obtiene 

d . 

xi X2— X2 xi = -t:{xi X2— X2 Xi) = , 

o tambien 

X'l X2— X2 xi = cte . 

Al sustituir t por f + T en la ecuacion anterior, el miembro derecho esta 
multiplicado por Hifi2 (debido a p8|); por lo que es evidente que 



/xi/i2 = 1 , (39) 

teniendo en cuenta (p7|) y sabiendo que los coeficientes son reales. Si x{t) es 
una integral de esta ecuacion, la funcion x * (t) tambien lo es. Lo anterior 
conduce a que /xi, H2 deben coincidir con nl, 112, es decir, o /ii=/^2 ^ tambien 
Ml y ^J'2 son reales. En el primer caso y teniendo en cuenta (|39|) resulta que 
Hi = 1/ fj-l lo que es igual |//i| = |At2| =1- En el segundo caso, las dos 
integrales son de la forma 

xi{t) = /u*/^F(t) 

X2{t) = /u"*/^G(t) . 

Una de estas funciones crece exponencialmente con el tiempo. 
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5. TRANSFORMACIONES CANONICAS 



Prologo: La idea principal de las transformaciones canonicas es encontrar 
todos aquellos sistemas de coordenas (en el espacio de fases) los cuales pre- 
serven la forma de las ecuaciones de Hamilton, independientem,ente de que 
Ham,iltoniano se trate. Posteriormente se escoge de entre todos esos sis- 
temas de coordenadas aquel en el cual se facilite la resolucion del problema 
en particular. 



CONTENIDO: 

5.1 Definiciones, Hamiltoniano y Kamiltoniano 

5.2 Condiciones necesarias y suficientes para que una transf. sea canonica 

5.3 Ejemplo de aplicacion de transf. canonica 
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5.1 Definiciones, Hamiltoniano y Kamiltoniano 

Se define una transformacion canonica, considerando los casos en que la 

transformacion depende explicitamente (o no) del tiempo, de la siguiente 

manera: 

Definicion 1: Una transformacion independiente del tiempo Q = Q{q,p), 

y P = P{q,p) se dice que es canonica si y solo si existe una funcion F{q,p) 

tal que 

dF{q,p) = ^Pidqi -^Pi{q,p)dQi{q,p) . 

i i 

Definicion 2: Una transformacion dependiente del tiempo Q = Q{q,p, t), y 
P = P{q,p, t) se dice que es canonica si y solo si existe una funcion F{q,p, t) 
tal que para un tiempo arbitrario fijo t = to 

dF{p,q,to) =^pidqi-^Pi{q,p,to)dQi{p,q,to) , 

i i 

donde 

dF{p, q, to) = 2^ g^^ dqi + }^ — dpi 

y 

.^, ,^ yr^dQ{p,q,to) , v-9Q(p,g,to) 

dQ{p, q, to) = 2^ ^ dqi + 2_^ dpi 

■ oqi . opi 

Ejemplo: Demostrar que la siguiente transformacion es canonica: 

Q = Tan~^ (^ 
\P 

Solucion: De acuerdo a la definicion 1, debemos verificar que pdq — PdQ 
es una diferencial exacta. Sustituyendo P y Q de la transformacion dada en 
el problema obtenemos: 

1^,2 , 2^Pdq-Qdq Jpq 



pdq - PdQ =pdq- -{p +q ) — :r- — ^ =d^ ^ 
2 p'^ + q'^ \ 2 

De lo anterior concluimos que efectivamente, la transformacion dada en el 
problema es canonica. Sabemos que un sistema dinamico se encuentra car- 
acterizado por su Hamiltoniano H = H{q,p,t), en donde q = q{qi, ^2, •••,<?«), 
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dH 


yi 


dpi 




dH 


Pi = 


dqi 



y p = pipi,P2, ■■■iPn), y que por tanto el sistema tiene asociado un conjunto 
de 2n ecuaciones diferenciales de primer orden, dadas por las ecuaciones de 
Hamilton: 

(1) 

(2) 

Sean las transformaciones de coordenadas en el espacio de fase, denotadas 
como 

Qj=Qj{q,P,t) (3) 

Pj = Pj{q,p,t) , (4) 

entonces, de acuerdo a lo dicho al principio, denotaremos como transforma- 
ciones canonicas al conjunto de transformaciones de la forma de (3) y (4) 
para las cuales, analogamente a (1) y (2), exista una funcion K = K{Q, P, t) 
tales que podamos escribir 

La relacion existente entre el Hamiltoniano H y el nuevo Kamiltoniano K^ 
se puede obtener a partir de las siguientes consideracionesQ 

De acuerdo al principio de Hamilton, la trayectoria real que un sistema 
clasico describira se puede obtener a partir de variacion de la integral de 
accion dada por 

6 j{Y,Pidq^-Hdt)=Q . (7) 

i 

Si la transformacion es canonica, entonces el nuevo Kamiltoniano K debe 
tambien cumplir una relacion similar a (7), es decir, con el nuevo conjunto 



'^Aqui seguiremos la misma terminologia empleada por Goldstein al referirse a la nueva 
funcion K = K{Q,P,t), la cual difiere del Hamiltoniano H = H(p,q,t) por una derivada 
temporal aditiva, como la Kamiltoniana. 

■^Una derivacion alternativa a la presente ha sido dada por G. S. S. Ludford and D. W. 
Yannitell, Am. J. Phys. 36, 231 (1968). 



75 



de coordenadas Q y P tambien es valido que 

5 j{Y,mQi-Kdt) = {). (8) 

i 

Sabemos ademas que, de acuerdo a la transformacion de Legendre, J2iPidqi— 
Hdt = L{q,q,t)dt, por lo que (7) -lo mismo que (8)- equivale a 

5 r L{q,q,t)dt = ^ (9) 

y que (9) no se altera si L es sustituido por C = L -\ ^|^ porque en este 

caso, 

5 t' Cdt = 5 f\L+'^^^^)dt, (10) 

Jti Jt± dt 

o en forma equivalente 

5 f ' Cdt = S f ' L{q,q,t)dt + 5Fiq(^2),t2) - SFiq(^,),ti) (11) 

, por lo que (10) y (11) difieren solamente por terminos constantes los cuales 
dan como resultado una variacion nula al moniento de aplicar el principio 
de Hamilton. 

De acuerdo a lo anterior podemos exigir que el Hamiltoniano H j el Kamil- 
toniano K se encuentren relacionados por la ecuacion |^ 

dF 
Piqi -H = P,Qi-K + — . (12) 

La funcion F es llamada funcion generadora. Esta puede ser expresada 
como una funcion de cualquier conjunto arbitrario de variables independi- 
entes. Sin embargo, algunos resultados muy convenientes son obtenidos si F 
es expresada como funcion de las n viejas variables y las n nuevas variables, 
mas el tiempo. Los resultados son especialmente convenientes si las n viejas 
variables son en su totalidad las n gj - o las n pj-, y si las nuevas variables 
son en su totalidad las n Qi - o las n Pj. 



^Algunos autores agregan al lado derecho de esta ecuacion un factor multiplicativo con- 
stante A ya que en este caso tampoco se ve alterado (1.9). Nosotros escogimos arbitraria- 
mente A = 1, es decir, decidimos trabajar con transformaciones canonicas restringidas, ya 
que este caso es suficiente para mostrar la estructura de las transformaciones canonicas. 
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Atendiendo a lo dicho en el parrafo anterior, las posibles combinaciones 
de n variables viejas y n variables nuevas -incluyendo a t- en la funcion 
generadora son: Q 

Fi = Fi{Q,q,t) (13) 

F2 = F2{P,q,t) 

F3 = Fs{Q,p,t) 

F4 = F4{P,p,t) . 

Por otro lado, si multiplicamos (1-12) por dt obtenemos: 

Pidqi - Hdt = PdQi - Kdt + dF . (14) 

Y haciendo el cambio F — > Fi en la relacion anterior, y recordando que dQi, 
dqi, y dt son variables independientes, obtenemos: 

dFi 



P 



K = H + 



dQi 
dFi 

dqi 

dFi 



dt 



Mediante manipuleos algebraicos es posible obtener expresiones analogas a la 
anterior que involucren a las restantes funciones generadores. Los resultados 
que se obtienen se muestran a continuacion: 

F4: Q.= It ,. = -11 K = B^^-§. 

En la practica resulta litil el siguiente teorenia, el cual permite, junto con 
las definiciones dadas en la introduccion para que una transformacion sea 
canonica (ya sea que la transforamacion dependa o no explicitaniente del 
tiempo), resolver cualquier problema mecanico de interes [|. 

^Usaremos la misma convencion de Goldstein para denotar cada una de las diferentes 
combinaciones de las variables nuevas y viejas en la funcion generadora. 
^Para un ejemplo, ver la seccion final de este capitulo. 

77 



Teorema 1.1.: Consideremos un sistema sobre el cual se ejerce una fuerza 
neta dada. Supongamos adenias que el estado dinamico del sistema esta 
definido por un conjunto de variables q,p = qi,q2, ■■■,Qn,Pi,P2, ■■■,Pn Y que 
el Haniiltoniano del sistema es H = H{q,p,t) tal que el comportamiento de 
las variables q y p este determinado por las ecuaciones de Hamilton 

dH{q,p,t) 



Pi 



dpi 
dH{q,p,t) 



dqi 
Si nosotros hacemos una transformacion a las nuevas variables 

Q = Q{q,p,t) ; P = P{q,p,t) 

y si la transformacion es canonica, es decir, si existe una funcion F{q,p,t) 
tal que para un tiempo arbitrario fijo t = to, 

dF{q,p, to) = Y^ Vidxi - ^ YidXi , 

i i 

donde Xi,yi = qi,pi o Pi,-qi y Xi,Yi = Qi,Pi, o Pi,-Qi, entonces las 
ecuaciones de movimiento en terminos de las variables Q y P son 

. _ dK{Q,P,t) 

^' ~ dP, 

dK{Q,P,t) 



p. = 

dQ 



donde 



j^_ jj , dF{q,p,t) ^^ dXi{q,p,t) 
K = H + + 2^ Yi- 



dt ^ ' dt 



Ademas si el determinante de la matriz [-q-^] es distinto de cero, entonces 
la ecuacion anterior se reduce a 

dt 

5.2 Condiciones necesarias y suficientes para que 
una transformacion sea canonica 

Hemos mencionado al principio que por transformacion canonica entender- 
emos aquella transformacion que, independientemente de cual sea la forma 
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del Hamiltoniano, preserva la forma de las ecuaciones de Hamilton. Se debe 
ser muy cuidadoso respecto a este punto ya que es posible que existan trans- 
fer maciones que preservan la forma de las ecuaciones de Hamilton, pero para 
un Hamiltoniano particular ^ Algunos autores le denominan a este tipo de 
transformacion transformacion canonica respecto a H [^ 

Para ilustrar este punto utilizaremos el siguiente ejemplo, sugerido en el 
articulo de J. Hurley: Consideremos un sistema fisico particular cuyo Hamil- 
toniano sea 

// = — (15) 

y consideremos las transformaciones 

Q : f. 

Es facil demostrar que el Kamiltoniano K dado por 

2p3/2 



K 



3m 



conduce a 



dK 

P = 2pp = = -—- 

oQ 

y 

. p pi/2 dK 
Q = q = — = = 7TF ■ ^^^^ 

m m oF 

Por otro lado, si escogemos el Hamiltoniano 

2m 

entonces es imposible encontrar un Kamiltoniano K tal que al usar las ecua- 
ciones de transformacion (16) se preserve la forma de las ecuaciones de 
Hamilton. Como vemos, en el anterior ejemplo las ecuaciones (16) preser- 
van la forma de las ecuaciones de Hamilton, pero para un Hamiltoniano en 
particular. 

®Ver por ejemplo, J. Hurley, Am. J. Phys. 40, 533 (1972). 

'^Ver por ejemplo, R. A. Matzner y L. C. Shepley, Classical Mechanics (Prentice Hall, 
1991). 
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Se puede demostrar que las condiciones necesarias y suficientes para que 
transformaciones de la forma (3) y (4) sean canonicas, es decir, que preser- 
ven la forma de las ecuaciones de Hamilton independientemente del Hamil- 
toniano que se considere, son: 

[Qi,Pj] = a (18) 

[n,Pj]=0 (19) 

[Qi,Qj]=0, (20) 

en donde a es una constante cualquiera, relacionada con cambios de es- 
cala. Algunos comentarios merecen ser mencionados antes de finalizar esta 
seccion. En primer lugar, debemos mantener en mente que Q y P no con- 
stituyen variables que definan la configuracion del sistema, es decir, no con- 
stituyen en general un conjunto de coordenadas generalizadas Q. Para dis- 
tinguir a Q y P, de las coordenadas generalizadas q y p, se les denomina 
variables canonicas. Y a las ecuaciones de movimiento -similares en forma 
a las ecuaciones de Hamilton para las coordenadas generalizas q y p- que 
se obtienen para Q y P se les denomina ecuaciones canonicas de Hamil- 
ton. En segundo lugar, aunque no lo probamos aqui, si la transformacion 
Q = QiQjPjt) y P = P^Q^Pt'I^) 6s canonica, entonces la transformacion in- 
versa q = q{Q, P,t) y p = p{Q, P, t) es tambien canonica ^ 

5.3 Ejemplo de aplicacion de TC 

Como se menciono en la introduccion a este capitulo, la idea principal de 

realizar una transformacion canonica es encontrar sistemas de coordenadas 

(en el espacio de fases) los cuales preserven la forma de las ecuaciones de 

Hamilton, independientemente de la forma del Hamiltoniano, y escoger de 

entre todas ellos, aquel que facilite la resolucion del problema en particular. 

Vamos a ilustrar este punto con un ejemplo. 

EJEMPLO : 

El Hamiltoniano de cierto sistema fisico esta dado por H = uj'^p{q + t)^, 

donde to es una constante. Determine q como una funcion del tiempo. 

Solucion: 



Salvo el caso trivial en que la transformacion canonica sea Q — q y P = p. 
^Para una demostracion sobre este punto, ver por ejemplo, E. A. Desloge, Classical 
Mechanics, Volume 2 (John Wiley & Sons, 1982) . 
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1. Resolviendo las ecuaciones de Hamilton para las variables q y p. Al 
aplicar las ecuaciones de Hamilton (1) y (2) al Hamiltoniano dado en este 
problema, se obtiene 

uj'^iq + tf=q, 2uj'^p{q + t) = -p . 

Este sistenia no es de solucion facil. Sin embargo, este problema se puede 
resolver facilmente con una adecuada transformacion canonica, tal como se 
muestra a continuacion. 

2. Haciendo uso de la transformacion canonica Q = q + t, P = p. De 
acuerdo al teorema dado en la seccion (1.1), puesto que 

op 

dP 

= 0, 



d{-c 
entonces el Kamiltoniano K del sistema esta dado por 

La forma de la funcion F{q, p, t) la encontramos a partir de la definicion de 
transformacion canonica dada en la seccion 1 -este caso corresponde a una 
transformacion canonica dependiente explicitamente del tiempo-. Es decir, 
a partir de 

dF{q,p,t) =pdq- PdQ . 

Y sustituyendo Q = q + t, P = p en la relacion anterior se obtiene sin 
ninguna dificultad que 

F{q,p,t) 

Por otro lado 



= c. 


c= constante . 


dP 

dt 


= 


dQ 

dt 


= 1 . 



Finalmente, sustituyendo estos resultados (y sustituyendo Q = q + t, P = p 
en el Hamiltoniano H) en (21) obtenemos 

K = P{uj^Q^ + 1) 
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Y de (5) obtenemos 



Q = u?Q^ + 1 . 



Esta ecuacion diferencial se resuelve facilmente, y se obtiene como resultado 

final 

1 

q = — tan(a;t + 0) — t , 

en donde ^ es una fase arbitraria. Como se puede observar, la correcta 
eleccion de la transformacion canonica puede facilitar la solucion de cualquier 
problema mecanico. 
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6. PARENTESIS DE POISSON 

Prologo: Los parentesis de Poisson son herramientas analiticas muy utiles 
para estudiar el comportamiento de cualquier sistenia dinaniico. Nosotros 
definiremos en este capitulo lo que se entiende por parentesis de Poisson; 
daremos algunas de sus propiedades, y finalmente presentaremos algunas 
aplicaciones de los mismos en el estudio de sistemas dinamicos. 



CONTENIDO: 

1. Definicion y propriedades 

2. Forniulacion de Poisson para las ecs. de movimiento 

3. Las constantes de movimiento en la forniulacion de Poisson 
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1. Definicion y propiedades de los parentesis de Poisson 

Si u y V son cualesquiera dos cantidades que dependen del estado dinamico 
del sistema (es decir, de p y de q) y posiblemente del tiempo, el parentesis 
de Poisson de u y v con respecto a un conjunto de variables canonicas q y p 
[^ es definido como 

r 1 ^r^ fdu{q,p,t)dv{q,p,t) du{q,p,t) dv{q,p,t)\ 

Los parentesis de Poisson tienen las siguientes propiedades (donde u, v, y w 
son funciones arbitrarias de q, p, y de t; a es una constante arbitraria, y r 
es cualquiera de las variables qi, pi o t) f^ 

Propiedad 1. [u-,v\ = — [^ji*] 

Propiedad 2. [u, u] = 

Propiedad 3. [u, v + w] = [u, v] + [u, w] 

Propiedad 4- [u,vw] = v[u,w] + [n, ujw 

Propiedad 5. a[u,v] = [au,v] = [u,av] 

Propiedad 6. « ^ [|^,^] + [^, |h] 

Propiedad 7. Identidad de Jacobi, [u, [v, w]] + [v, [w, u]] + [w, [u, v]] = . 

Otra propiedad muy importante de los parentesis de Poisson es enunciada 
en el siguiente teorema: 

Teorema 6.1 Si la transformacion Q = Q{q,p,t), P = P{q,p,t) es una 
transformacion canonica, el parentesis de Poisson de dos cantidades u y v 
con respecto al conjunto de variables q, p, es igual al parentesis de Poisson 
de u y V con respecto al conjunto de variables Q, P, es decir 

y^ f du{q,p,t) dv{q,p,t) _ du{q,p,t) dv{q,p,t) 
'Y \ dqi dpi dpi dqi 



^"Como en el capitulo anterior, entendemos a q y p como q — gi,g2,---,qn y p = 
pi,P2,...,Pn. 

'^^La prueba de estas propiedades se logra utilizando la definicion de los parentesis de 
Poisson para expresar cada uno de los terminos de estas identidades en terminos de las 
derivadas parciales de u, v, y w, j notando por inspeccion la valides de las ecuaciones 
resultantes. 



v^ fdu{q,p,t) dv{q,p,t) du{q,p,t) dv{q,p,t) 



V dQi dP^ dPi dQ^ 



2. La formulacion de Poisson para las ecuaciones de movimiento 

Resumimos a continuacion, en forma de teoremas, los resultados mas im- 
portantes sobre los parentesis de Poisson en el analisis de el movimiento de 
cualquier sistema dinamico P^ : 

Teorema 6.2: Consideremos un sistema cuyo estado dinamico esta definido 
por las variables canonicas q^py cuyo comportamiento dinamico esta definido 
por la Hamiltoniana H = H{q,p,t). Sea F una cantidad arbitraria que de- 
pende del estado dinamico del sistema (es decir, de q, p, y posiblemente t). 
La razon de cambio en el tiempo de F esta dada por 



dt 

donde [F, H] es el parentesis de Poisson de F con H. 

Teorema 6.3 (Formulacion de Poisson de las ecs. de movimiento). 

Consideremos un sistema cuyo estado dinamico esta definido por las vari- 
ables canonicas q, p, y cuyo comportamiento dinamico esta dado por la 
Hamiltoniana H = H{q,p,t). El movimiento del sistema esta gobernado 
por las ecuaciones 

4i = [Qi,H] 
Pi = [Pi,H] . 

3. Las constantes de movimiento en la formulacion de Poisson 

Nuevamente enunciaremos los resultados mas importantes sobre la formu- 
lacion de Poisson de las constantes de movimiento, en forma de teoremas. 
Estos son los siguientes: 

Teorema 6.4 Si una cantidad dinamica F no es una funcion explicita de el 
tiempo, y si el parentesis de Poisson de F con H es nulo, es decir, [F, H] = 0, 
entonces F es una constante del movimiento. 

Corolario 6. 4. a Si el Hamiltoniano no es una funcion explicita del tiempo, 
entonces es una constante de movimiento. 



^^ Las demostraciones han sido omitidas por ser muy conocidas. Ver por ejemplo, E. 
A. Desloge, Classical Mechanics, Volume 2 (John Wiley & Sons, 1982). 
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7. LAS ECUACIONES DE HAMILTON- JACOBI 



Prologo: Sabemos de los capitulos anteriores que podemos, en principio, 
disminuir la complejidad de cualquier problema de dinaniica escogiendo una 
adecuada transformacion canonica. En particular, podemos tratar de buscar 
aquellas transforniaciones canonicas que hagan que el Kamiltoniano K sea 
nulo, lo que da las ecs. de Hamilton-Jacobi. 



CONTENIDO: 

7.1 Intro duccion 

7.2 Ec. de Hamilton-Jacobi dependiente del tiempo 

7.3 Ec. de Hamilton-Jacobi independiente del tiempo 

7.4 Generalizacion de la ec. de Hamilton-Jacobi 

7.5 Ejemplo de aplicacion de la ec. de Hamilton-Jacobi 
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7.1 Introduccion. 

Para lograr nuestros propositos en este capitulo, necesitamos hacer uso del 
siguiente resultado, el cual nos permite encontrar el conjunto de variables 
canonicas que hacen que el Kaniiltoniano asuma una forma particular. 
Teorema 7.1. Consideremos un sistema cuyo estado dinaniico esta definido 
por las variables p, q y cuyo comportamiento bajo la accion de una fuerza 
dada esta dado por la Hamiltoniana H = H{q,p,t). Sea K = K{Q,P,t) 
una funcion conocida de las variables canonicas Q, P, y tambien del tiempo. 
Entonces cualquier funcion F{q, Q, t) que satisfaga la ecuacion diferencial 
parcial 



K 



dFiq,Q,t) ■ 
.^' dQ '\ 



H 



dF{q,Q,t) 



y tambien satisface la condicion 

d^F{q,Q,t) 



dqjdQj 



dq 



/O 



,t 



+ 



dF{q,Q,t) 
dt 



es una funcion generadora para una transformacion canonica de las variables 
q, p a las variables Q, P, y la correspondiente Kamiltoniana es justamente 
K = K{Q,P,t). 

En las siguientes secciones haremos uso del anterior problema para encontrar 
aquellas transformaciones canonicas cuyo kaniiltoniano sea nulo |^; esto nos 
conducira a las ecuaciones de Hamilton-Jacobi. 



7.2 Ecuacion de HJ dependiente del tiempo. 

Como una consecuencia del Teorema 7.1 y de exigir que el Kaniiltoniano K 
sea nulo, obtenemos el siguiente teorema: 

Teorema 7.2. Consideremos un sistema de / grados de libertad cuyo estado 
dinaniico este definido por el conjunto de variables q, py cuyo Hamiltoniano 
sea H = H{q,p,t). Si construimos la ecuacion diferencial parcial 



H 



dSiq,t) 
dq 



+ 



dSiq,t) 
dt 







(1) 



^Es decir, K [Q, 



BQ ' J ^ "■ 
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y si encontrartios una solucion a esta ecuacion de la forma 

S = S{q,a,t) , 

donde a = ai, a2, ■■■, «/ es un conjunto de constantes, y si la solucion satis- 
face la condicion 

d'^S{q,a,t) 



/O 



dqidai 

entonces q{t) puede ser obtenido de las ecuaciones 

dS{q,a,t) 



dui 



Pi, 



donde /3 = (3i,(32, ■■■,f3f es un conjunto de constantes. El conjunto de ecua- 
ciones (2) nos proporciona / ecuaciones algebraicas en las / variables de- 
sconocidas qi,q2, ■■■,qf. Los valores de las constantes ay (3 son determinados 
por las condiciones de frontera. Ademas es posible, dada q{t), encontrar p(t) 
a partir de la relacion 



Pi 



dS{q, a,t) 
dqi 



(3) 



A la ecuacion diferencial parcial (1) se le denomina ecuacion de Hamilton- 

Jacobi dependiente del tiempo. Y a la funcion S{q, a, t) se le denomina 

funcion principal de Hamilton. 

Para lograr una mayor comprension de este teorema, asi como del signifi- 

cado de las constantes a y (3, procedemos a continuacion a desarrollar su 

demostracion. 

Demostracion del Teorema 7.2. De acuerdo al Teorema 7.1, cualquier 

funcion F(q, Q, t) que satisface la ecuacion diferencial parcial 



H 



dF{q,Q,t) 
dq 



dF{q,Q,t) ^^ 



dt 



y tambien satisface la condicion 

d^F{q,Q,t) 



dqjdQj 



/O 



debera ser una funcion generadora para una transformacion canonica a un 
conjunto de variables canononicas Q, P para las cuales el Kamiltoniano K 
sea nulo, es decir K{Q, P, t) = 0. La funcion 



F{q,Q,t) = [S{q,a,t)]a=Q = S{q,Q,t) 



es tal funcion. Entonces S{q, Q, t) es la funcion generadora para una trans- 
fer macion canonica que lleva al nuevo conjunto de variables Q, P, para las 
cuales el Kaniiltoniano K es identicamente igual a cero. Las ecuaciones de 
transformacion asociadas con la funcion generadora S{q, Q, t) estan definidas 
por las ecuaciones 

dS{q,Q,t) 

Pi = ^ (4) 

dqi 
_ dSiq,Q,t) 

^' ~ dQ— ^^^ 



dK(Q,P,t) 
^P^ 



y conio K(Q, P, t) = 0, las ecuaciones de movimiento son 

= 

p _ dK{Q,P,t) _ p. 

De las anteriores ecuaciones inferimos que 

Qi = OLi (6) 

P^ = -/3, (7) 

donde «» y /?« son constantes. La eleccion el signo negativo para /? en 
la ecuacion (7) se hizo por mera conveniencia. Si ahora sustituimos las 
ecuaciones (6) y (7) en la ecuacion (5) obtenemos 



-A 



dS{q,Q,t) 
dQi 



dS{q, a,t) 
Q=a dui 



la cual se reduce a (2). Si ademas sustituimos (6) en (4) obtenemos (2). 
Esto completa la prueba. 

7.3 Ecuacion HJ independiente del tiempo. 

Si la funcion Hamiltoniana no depende explicitamente del tiempo, nosotros 
podemos parcialmente resolver la ecuacion de Hamilton-Jacobi dependiente 
del tiempo. Este resultado se enuncia en el siguiente teorema: 
Teorema 7.3. Consideremos un sistema de / grados de libertad cuyo es- 
tado dinamico esta definido por el conjunto de variables q, p, y cuyo compor- 
tamiento bajo la accion de una fuerza definido por la funcion Hamiltoniana 
independiente del tiempo H{q,p). 

89 



Si construimos la ecuacion diferencial parcial 

dW{q) 



H 



dq 



E, (8) 



donde E es una constante cuyo valor para un conjunto particular de condi- 
ciones es igual al valor de la constante de movimiento H{q,p) para las condi- 
ciones de frontera dadas, y si nosotros podemos encontrar una solucion a 
esta ecuacion, de la forma 

W = W{q,a) , 

donde a = ai,a2, ■■■, a/ es un conjunto de constantes que explicita o implicitamente 
incluye a la constante E, es decir, E = E{a), y si la solucion satisface la 
condicion 

d'^Wiq,a) 



/O, 



dqidaj 
entonces las ecuaciones de movimiento estan dadas por 

dS{q,a,t) 



donde 



A (9) 



S{q, a, t) = W{q, a) - E{a)t 

y (3 = Pi, (32, ■■■,(3f es un conjunto de constantes. El conjunto de ecuaciones 
(9) proporciona / ecuaciones algebraicas en las / variables desconocidas 
qi,q2, ■■■,qf. Los valores de las constantes a y (3 son determinados por las 
condiciones de frontera. La ecuacion diferencial parcial (8) recibe el nombre 
de ecuacion de Hamilton- J acohi independiente del tiempo, y a la funcion 
W{q, a) se le conoce como funcion caracteristica de Hamilton. 

7.4 Generalizacion de la ecuacion de HJ. 

La ecuacion de Hamilton-Jacobi puede ser generalizada como se establece 
en el siguiente teorema, el cual permite, algunas veces, el simplificar algunos 
problemas de Hamilton-Jacobi. 

Teorema 7.4. Consideremos un sistema de / grados de libertad cuyo estado 
dinamico esta definido por un conjunto de varibles x, y, donde Xi,yi = qi,Pi 
o Pi, —qi, y cuyo comportamiento bajo la accion de una fuerza dada esta 
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dado por la Hamiltoniana H{x,y,t). Si construimos la ecuacion diferencial 
parcial 



H 



dS{x,t) , 
ox 



, dS{x,t) _^ 
dt 



y si podemos encontrar una solucion a esta ecuacion de la forma 

S = S{x, a,t) 

donde a = ai, a2, ■■■, cuf es un conjunto de constantes, y si la solucion satis- 
face la condicion 

d'^S(x,a,t) 



dxjdaj 



/O 



entonces el movimiento del sistema puede ser obtenido de las ecuaciones 

dS{x,a,t) 



dxi 
dS{x,a,t) 



Pi 



(10) 
(11) 



donde (3 = Pi, (32, ■■■, Pf es un conjunto de constantes. 

En la siguiente seccion mostrarenios algunos ejemplos donde hacemos uso 

de todos los resultados obtenidos hasta aqui en este capitulo. 



7.5 Ejemplo de aplicacion de la ecuacion de HJ 

Vamos a resolver el problenia del oscilador armonico unidimensional de masa 
m, utilizando el metodo de Hamilton-Jacobi. 
Sabemos que la Hamiltoniana del sistema es 



^^ jP_ kx^ 
2m 2 



(12) 



De acuerdo al Teorema 7.2 la ecuacion de Hamilton-Jacobi para el sistema 

es 

1 /dF\ kq^ dF , , 



2m\dq)^ 2 ^ dt ° 



Supongamos una solucion a (13) de la forma F = Fi[q) + F2{t). Por lo 
tanto, (13) se convierte en 



2m\dq ) 2 



dF2 
dt 



(14) 
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Haciendo cada lado de la ecuacion anterior igual a a, encontramos 



1 (dF,y kg- _ ^ 
2m\dq J ' 2 


(15) 


dF2 

dt ~ "" 


(16) 


Omitiendo las constantes de integracion, las soluciones son 




Fi = J^2mia-^)dq 


(17) 


Fa = -at 


(18) 



Por tanto, la funcion generadora F es 



F = / J2m(a - —)dq - at (19) 



Y de acuerdo a (2), q(t) se obtiene a partir de 



< / J2m{a - -^)dq - at > 



(3 = ^{ j \l2m[a--^)dq-at) (20) 



2m f dq 



a-T 



Y realizando la integral obtenemos que 



Y despejando q obtenemos finalmente que 



(21) 



'^sen''\qJk/2a) = t + (3 (22) 

k * 



sen^k/m{t + /?) . (23) 



k 

Podemos ademas dar una interpretacion fisica a la constante a siguiendo el 
siguiente razonamiento: 

El factor y-^ debe corresponder justamente a la maxima amplitud A que 
el oscilador puede tener. Por otro lado, la energia total E de un oscilador 
armonico undimensional de amplitud A esta dada por 

^ 1, ,2 1, / f^V 
E = -kA^ = -k\\ — =a. 
2 2 \\l k 
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Es decir, a es fisicamente igual a la energia total E del oscilador armonico 
unidimensional. 
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8. VARIABLES ACCION-ANGULO 



Prologo: La ecuacion de Hamilton-Jacobi nos proporciona un metodo de 
transformar un conjunto de variables q, p aun segundo conjunto de variables 
canonicas Q, P, siendo cada una de las variables canonicas una constante 
de movimiento. 

Ahora, en este capitulo nosotros consideraremos un metodo, valido para 
ciertos tipos de movimiento, para transformar un conjunto de variables q, 
p, a un segundo conjunto de variables Q, P, no siendo al mismo tiempo las 
dos variables canonicas constantes de movimiento. 



CONTENIDO: 

8.1 Sistemas separables 

8.2 Sistemas ciclicos 

8.3 Variables accion-angulo 

8.4 Movimiento en variables accion-angulo 

8.5 Importancia de las variables accion-angulo 

8.6 Ejemplo: el oscilador armonico 
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8.1 Sistemas separables. 

Por sistemas separables nosotros entenderemos sistemas cuya Hamiltoniana 
no sea una funcion explicita del tiempo, es decir 

H = H{q,p) 

y para los cuales es posible encontrar una solucion a la ecuacion de Hamilton- 
Jacobi independiente del tiempo, de la forma 

W{q,a)=J2W^{q^,a) . 



8.2 Sistemas ciclicos. 

Sabemos que el estado de un sistema se encuentra caracterizado por un 
conjunto de coordenadas generalizadas q = qi,q2, ■■■iQf Y por el conjunto 
de momenta generalizados p = pi,p2, ■■.,Pf- Conforme el sistema se mueve, 
este traza una orbita en el espacio q, p (el espacio cuyas coordenadas son 
las variables qi,q2, ■.■,qf,pi,P2, ■■-,?/)■ De igual forma, el sistema traza una 
orbita en cada uno de los subespacios qi,Pi- La orbita en cada uno de los 
pianos qi,Pi puede ser representada por una ecuacion de la forma pi = Pi{qi) 
o por un par de ecuaciones pi = Pi{t), qi = qi(t). Si para cada valor de 
i, la orbita pi = Pi{qi) en el piano qi — pi es una curva cerrada, entonces 
nos referimos a tal sistema como un sistema ciclico. En la siguiente figura 
mostramos las dos posibilidades por las cuales un sistema puede ser ciclico. 
En la Figura 8.1a, el sistema es ciclico porque qi oscila entre los limites 
definidos qi = a y qi = b; en la Figura 8.16, el sistema es ciclico porque qi se 
mueve desde qi = a hasta qi = b, y luego inicia nuevamente en qi = a. 
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Es necesario en este punto hacer dos importantes aclaraciones: 

Aclaracion 1: El termino ciclico ha sido introducido solamente para 
simplificar la notacion en las siguientes secciones. No debe interpretarse este 
termino pensando que si el sistema es ciclico en cada uno de los subespacios 
qi,pi, entonces el sistema deba regresar a su estado original en el espacio 
q,p. 

Aclaracion 2: Si un sistema ciclico solamente tiene un grado de liber- 
tad, el tiempo requerido por el sistema para ejecutar un ciclo en el piano 
q — p es constante; por tanto el movimiento en el piano q — p sera periodico 
en el tiempo. Si el sistema tiene mas de un grado de libertad, entonces, 
en general, el tiempo requerido para ejecutar un ciclo particular en uno de 
los espacios qi,Pi no sera una constante, sino que dependera del movimiento 
de las otras coordenadas; por tanto, el movimiento en el espacio dado qi,Pi 
no sera periodico en el tiempo. Debe de tenerse cuidado en este punto, ya 
que no es cierto que cada uno de los movimientos en cada piano qi,Pi es 
periodico en el tiempo. 

8.3 Variables accion-angulo. 

Consideremos un sistema ciclico de / grados de libertad, cuyo estado dinamico 
esta caracterizado por las variables canonicas q, p. Sea H{q,p) la Hamilto- 
niana del sistema, y sea 



W{q,a)^Y.^i(li^ 



a] 



(donde a = ai,a2, ■■■,af son constantes) una solucion a la ecuacion de 
Hamilton-Jacobi independiente del tiempo 

dW 

Sea J = Ji, J2, ..., J/ el conjunto de constantes definidas por las ecuaciones 

Ji{a) = f ^- dqi , (1) 

donde la integral es sobre un cilo completo para la variable qi. Si nosotros 
usamos la funcion 

W{q,a) = W[q,a{J)] 
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Pi 



i 

^ Y^wMuJ) 

i 

como una funcion generadora para una transformacion canonica de las vari- 
ables g, p a un nuevo conjunto de coordenadas w, y momenta J, es decir, si 
definimos las variables w j J por las ecuaciones de transformacion 

dW{q,a) ^ dWi{qi,J) 
dqi dqi 

dW{q,J) 

entonces las nuevas coordenadas wi,W2, ■■■,Wf son llamadas variables dngulo, 
y los nuevos momenta Ji, J2, ..., J/ son llamadas variables accion. 
De la ecuacion (2) obtenemos que 

P.te,a) = ^^%^ = ^^^^. (4) 

dqi dqi 

Sustituyendo la ecuacion (4) en (1) obtenemos 

Ji{a) = f pi{qi,a)dqi . (5) 



La ecuacion pi = pi{qi, a) es la ecuacion de la proyeccion de la orbitap = p{q) 
sobre el subespacio pi,qi. La integral en el lado derecho de la ecuacion (5) 
es por tanto el area encerrada dentro de la orbita, o bajo la orbita, como se 
ilustra por las regiones sombreadas en la Figura 8.1. Por tanto, la funcion 
Ji{a) puede ser interpretada geometricamente como el area barrida en el 
subespacio qi,Pi durante un ciclo completo en este subespacio. Esta area 
depende de las constantes a o equivalentemente de las condiciones iniciales 
y puede en general asumir cualquier valor pi. 

8.4 El movimiento de un sistema en terminos de las 
variables de accion-angulo. 

Resumimos el movimiento de un sistema en terminos de las variables de 
accion-angulo en la siguiente teorema: 

^^Historicamente, el primer intento de pasar de la mecanica clasica a la mecanica 
cuantica consistio en asumir que el valor que Ji podia tomar no era completamente arbi- 
trario, sino que debe ser un miiltiplo de h/2n, donde h es la constante de Planck. 
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Teorema 8.4. 

Considerenios un sistenia ciclico separable de / grades de libertad cuyo 
estado dinamico esta dado por las variables q,p = qi,q2, ■■■,qf,Pi,P2, ■■-,?/ , 
y cuyo comportamiento dinamico esta dado por el Hamiltoniano H{q,p). Si 
nosotros transformamos a las variables de accion-angulo J, w, entonces el 
Hamiltoniano H es una funcion de J solamente, es decir, 

H = H{J) 

y el movimiento del sistema esta dado por 

Ji = li 
Wi = Vit + 4>i , 

donde 7^ y (pi son constantes que estan determinadas por las condiciones 
iniciales, y las Vi son constantes, llamadas las frecuencias del sistema y se 
encuentran definidas como 



I'i 



dH{J) 



J=ii 



8.5 Importancia de las variables accion-angulo 

La importancia de las variables accion-angulo radica en que proporcionan 
una tecnica potente para la obtencion de la frecuencia de un movimiento 
periodico de un movimiento sin hallar una solucion completa del movimiento 
del sistema. 

Lo anterior puede ser visto del siguiente argumento: Consideremos el cambio 
de w cuando q describe un ciclo completo, dado por 

I dw , 
i\w = <p -T^dq 
J oq 

Por otro lado, sabemos que 

dW 

por lo que 



d fdW^ 
dq 



dJ J dq 
d r 

= 1 . 

El anterior resultado nos indica que w cambia en una unidad cuando q varia 
a lo largo de un periodo completo. 
De la relacion 

w = ft + 4> , 

concluimos que en un periodo r 

Alt; = 1 

= VT , 

es decir, podemos identificar a la constante v con el reciproco del periodo, 

1 

u = — . 

T 

8.6 Ejemplo: El oscilador armonico simple 

A partir del formalisto de las variables accion angulo vamos a demostrar que 
la frecuencia u del oscilador armonico simple unidimensional esta dada por 

Dado que H es una constante de movimiento, la orbita en el espacio q — p 
esta dada por 

p^ kq^ ^ 
2m 2 



donde E es la energia. Esta es la ecuacion de una elipse con semiejes \/2mE 
y \j2Elk. El area encerrada por la elipse es igual al valor de la variable de 
accion J. Por tanto, 

I \2E [m 

V K \ k 

Se sigue entonces que 



H(J) = E = ySp.j 
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y la frecuencia esta dada por 



dH{J) _ ^fkj^ 
dJ ~ 27r 
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9. TEORIA CANONICA DE PERTURBACIONES 



Prologo: Existen muchos problemas en la naturaleza que no pueden ser 
resueltos de manera exacta. Por esta razon, y tomando en cuenta el gran 
desarrollo experimentado por la informatica, se ha puesto mucho interes en 
el desarrollo de metodos para hallar soluciones aproximadas. El metodo de 
perturbaciones se aplica cuando se tiene un problema fisico que no se puede 
resolver exactamente, pero cuya Hamiltoniana difiere solo ligeramente de 
la Hamiltoniana correspondiente a un problema que puede resolverse de 
manera exacta. A la diferencia entre ambas Hamiltonianas se le conoce 
como la Hamiltoniana de la perturbacion y la teoria de perturbaciones esta 
basada en la pequiiez de la misma. 



CONTENIDO: 

9.1 Teoria de perturbaciones dependiente del tiempo (con dos ejemplos) 

9.2 Teoria de perturbaciones independiente del tiempo (con un ejemplo) 
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9.1 Teoria de pertubaciones dependiente del tiempo 

La formulacion de la mecanica clasica que simplifica mas el desarrollo de la 
teoria de perturbaciones es la de Hamilton-Jacobi. Consideremos entonces 
que Hf)(p,q,t) es la Hamiltoniana correspondiente al problema soluble o no 
perturbado, y que se ha solucionado mediante la funcion principal de Hamil- 
ton S{q,ao,t), la cual genera una transformacion canonica, de las coorde- 
nadas {p,q) a (aoj/^o)) sn la que la nueva Hamiltoniana (o Kamiltoniana) 
Ko del sistema no perturbado, es nula. En simbolos: 



— + Ho{ — ,q,t)=Ko = 0. (1) 

Esta es la ecuacion de Hamilton-Jacobi y hemos usado que p = dS/dq. Las 
coordenadas canonicas transformadas (ao, /3o) son entonces todas constantes 
en el caso no perturbado, ya que Kq = y: 






dKo 

"5/3o' 
dKo 



dao 
Consideremos ahora la hamiltoniana del sistema perturbado como: 

Hiq,p,t)=Ho{q,p,t)+AH{q,p,t); (AH ^ Hq). (3) 



Aunque (ao,/3o) siguen siendo coordenadas canonicas (pues la transfor- 
macion generada por S es independiente de la forma particular de la hamil- 
toniana), ya no seran constantes y la Kamiltoniana del sistema perturbado 
(K) no sera nula. Para no olvidar que en el sistema perturbado las coor- 
denadas transformadas ya no son constantes, las denotaremos como a y /3, 
en vez de ao y /3o) que seri'an sus correspondientes valores constantes en el 
sistema no perturbado. En el sistema perturbado tenemos pues que: 



Kia, f3,t) = H + — = iHo + —) + AH = AH{a, 13, t). (4) 



Las ecuaciones de movimiento que satisfacen las variables transformadas del 
sistema pertubado son entonces: 
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a,: 



A 



dAH{a,p,t) 
dAH{a,f3,t) 



(5) 



donde i = 1,2, ...,n y n es el numero de grades de libertad del sistema. 
Las ecuaciones anteriores son rigurosas. Si del sistema de 2n ecuaciones se 
pudieran obtener ai y (3i como funcion del tiempo, las ecuaciones de trans- 
formacion {p, q) -^ (a, (3) darian pi y qi en funcion del tiempo y el problema 
estaria resulto. Sin embargo, la solucion exacta de las ecuaciones en (5) 
no suele ser menos dificil que la solucion de las ecuaciones originales. De 
(5) vemos que aiin cuando a y /3 ya no son constantes, su variacion con el 
tiempo es lenta, si suponemos que /S.H cambia poco respecto a a y (3. Una 
primera aproximacion a la variacion temporal de (a, f3) se obtiene susti- 
tuyendo en los segundos miembros de (5) a a y (3 por sus valores constantes 
o no perturbados; es decir: 



an 
Pii 



dAH{a,(3,t) 



d(3i 
dAH{a,(3,t) 



dai 



(6) 



donde an y Pn representan las soluciones a primer orden de perturbacion 
para ai y (3i, y las barras verticales con subi'ndece cero indican que despues 
de la derivacion deben sustituirse a y (3 por sus valores constantes no per- 
turbados. Hecho esto, las ecuaciones en (6) se pueden integrar para dar las 
ai y las /3j en funcion del tiempo (a primer orden). Luego, mendiante las 
ecuaciones de transformacion se obtienen p y q como funcion del tiempo en 
una primera aproximacion. La aproximacion de segundo orden se obtiene 
sustituyendo en los segundos miembros de (6) la primera aproximacion de 
la dependencia de a y f3 con respecto al tiempo. En general, la solucion de 
perturbacion de orden N se obtiene integrando las ecuaciones: 





djAT = 


dAH{a,p,t) 
df3i 




PiN = 


dAH{a,p,t) 
dai 


Ejemplo 1 







N-l 



N-1 



(7) 



Consideraremos aqui el caso simple de una particula libre, que mas adelante 
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sujetaremos a una perturbacion armonica simple. Este ejemplo, aunque triv- 
ial, servira para ilustrar el procedimiento delineado antes. La Hamiltoniana 
no pertubada es: 






Puesto que Hq ^ Hq{x), es decir, como x es ciclica, p = a^ es una constante 
en el sistema no perturbado. Recordando que p = dS/dx y sustituyendo en 
(1): 



'^V + - = (9) 

2m \dx) ^ dt ^ ^ 



Ahora bien, dado que el sistema es conservativo, resulta conveniente consid- 
erar una funcion principal de la forma: 



S = Six)+Fit). (10) 



Esta clase de separacion de variables, es especialmente litil cuando la Hamil- 
toniana no depende explicitamente del tiempo, en donde se propone que 
F{t) = —Et, con E como la energia total del sistema ^ Con (10) en (9) 
tenemos que: 

1 / y-7G \ ^ 

(--) =E, S = \/2mEx = aox. (11) 

2m \dxj 

Y sustituyendo (11) en (10), junto con el hecho de que en este caso la 
Hamiltoniana es igual a la energia, tenemos que la funcion principal de 
Hamilton es: 



S = aox--—. (12) 

Im 



Si el momentum transformado es oq, la coordenada transformada (que 
tambien es constante en el sistema no perturbado) es: 



dS aot 

Po = t; — = x 

oan m 



'Vease: Spiegel, Murrary R. Mecdnica Teorica, pp. 315, 316. 
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de modo que la transformacion generada por S esta dada por las ecuaciones: 



P = "0, 



^ + /?o, (13) 

m 



que es la solucion esperada para el movimiento de una particula libre. Lo 
realizado hasta aqui solo muestra el procedimiento para hallar las ecuaciones 
del movimiento mediante la formulacion de Hamilton- Jacobi. Es hasta ahora 
que introduciremos una perturbacion de la forma: 



AH=^ = ^^, (14) 



o bien, en terminos de las coordenadas transformadas, utilizando (13): 



mu;^ /'at \^ , , 

AH = ^-(- + p) . 15 

2 \m J 



Notese que en la expresion anterior hemos suprimido ya los submdices de 
las coordenadas transformadas, pues estamos ya ocupandonos del sistema 
perturbado. 

Sustituyendo (15) en (5): 



2 /«^ 
a = —muj \- p 

p = u;h(- + p). (16) 



Las ecuaciones anteriores tienen una solucion exacta y es de forma armonica, 
como cabe esperar. Para asegurarnos solo derivamos respecto al tiempo la 
primera de las ecuaciones y Uegaremos a que a tiene una variacion armonica 
siemple, lo que tambien puede asegurarse de x en virtud de las ecuaciones 
de transformacion dadas en (13), que siguen manteniendo su forma en el sis- 
tema perturbado (omitiendo, claro, los submdices de las coordenadas trans- 
formadas). Sin embargo, nos interesa ilustrar el metodo de perturbaciones, 
asi que consideremos que k (la constante elastica) es un parametro pequeiio 
y busquemos soluciones aproximadas de distintos ordenes de perturbacion, 
sin perder de vista que las variables transformadas (a, (3) en el sistema per- 
turbado dejaran de ser constantes. Dicho de otra manera, aiin cuando (a, (3) 
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contienen informacion referente a los parametros del sistema no perturbado, 
el efecto de la perturbacion es hacer variar estos parametros con el tiempo. 

La perturbacion de primer orden se obtiene segiin esta indicado de manera 
general en (6). Asi que debemos sustituir en los segundos miembros de (16), 
a y (3 por sus valores no perturbados. Para simplificar consideremos que 
x{t = 0) = y por tanto que (3q = 0, entonces: 



Oil = —^ Ciot, 
■n ^'^^^ 

Pi = ao , (17) 

m 



que luego de integrar nos conduce a: 



«! = ao 



w^aot^ 



A = ^^ (18) 

3m 



Las soluciones para x y p a primer orden las obtenemos sustituyendo ai y 
Pi en las ecuaciones de transformacion (13), de donde: 



^ (cot - ^^ 



mto 



coH^ 



p = ao 1 - ^— . (19) 



Para generar la solucion aproximada a un segundo orden de perturbacion 
debemos hallar 0:2 y /32 5 segiin se indico en (7), sustituyendo en los segundos 
miembros de (16), a y (3 por ai y (3i como fueron dadas en (18). Integrando 
CK2 y /?2 y luego utilizando de nuevo las ecuaciones de transformacion (13), 
llegamos a las soluciones de segundo orden para x y p: 



muj \ 3! ^' / 

P = «o 1-^ + ^ . (20) 



En el limite en que el orden de perturbacion N tiende a infinito, obtenemos 
las soluciones esperadas y compatibles con las condiciones iniciales: 
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x^ sinwt, p^aocosojt. (21) 

Las variables transformadas (q, (3) contienen informacion referente a los 
parametros de la orbita sin perturbar. Por ejemplo, si consideramos como 
sistema no perturbado aquel correspondiente al problema de Kepler, un sis- 
tema adecuado de coordnadas (a, (3) podrian ser las variables (J, 5) que son 
respectivamente la variable accion y el angulo de fase que aparece en en 
la variable angulo w (recordar que w = vt + 5, donde v es la frecuencia). 
Estas variables estan relacionadas con los parametros orbitales tales como 
el semieje mayor, la excentricidad, la inclinacion, etc. 

El efecto de perturbacion es hacer variar estos parametros con el tiempo. Si 
la perturbacion es pequeiia, la variacion de los parametros durante un periodo 
del movimiento no perturbado tambien sera pequeiia. De modo que por 
cortos intervalos de tiempo el sistema se movera a lo largo de una orbita, 
llamada orbita osculatriz, que es de la misma forma funcional que la del 
sistema no perturbado; sin embargo, los parametros de esta orbita vari'an 
con el tiempo. 

Los parametros de la orbita osculatriz pueden variar con el tiempo de dos 
maneras: 

• Variacion periodica: el parametro vuelve a su valor inicial despues 
de un intervalo de tiempo que en primera aproximacion suele ser el 
periodo del movimiento no perturbado. Estos efectos de la pertur- 
bacion no alteran los valores medios de los parametros y por tanto la 
trayectoria sigue siendo muy parecida a la orbita no perturbada. Es- 
tos efectos se pueden eliminar promediando la perturbacion sobre un 
periodo del movimiento no perturbado. 



Variacion secular: Al final de cada uno de los periodos orbitales suce- 
sivos hay un incremento neto del valor del parametro. Al cabo de 
muchos periodos los parametros orbitales pueden ser muy diferentes 
de sus valores no perturbados. Rara vez interesa el valor instantaneo 
de la variacion de algiin parametro, digamos la frecuencia, porque su 
variacion es muy pequeiia en la mayoria de casos en que funciona 
el formalismo de la teoria de las perturbaciones. (Esta variacion es 
tan pequeiia que resulta dificil, sino imposible, percibirla en uno solo 
periodo orbital, y por eso se mide solo la variacion secular despues de 
varies periodos.) 



Ejemplo 2 
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Del problema de dos cuerpos tenemos que si al potencial de Kepler se le 
suma un potencial de la forma 1/r^, la orbita del problema acotado es una 
elipse que gira y cuyo periapsis esta animado de precesion. Encontraremos 
aqui la velocidad de precesion, considerando un potencial perturbador algo 
mas general: 



1- = ----, (22) 



donde n(> 2) es un entero, j hes tal que el segundo termino del potencial sea 
una pequeiia perturbacion del priniero. La Haniiltoniana de la perturbacion 
sera pues: 



^H = -^. (23) 



En el problema sin perturbar, la posicion angular del periapsis en el piano 
de la orbita viene dada por la constante to = 2itw2- En el caso perturbado: 



. ^ dAH dAH ,^^, 

donde hemos usado J2 = 2ttI. Ademas J2 y W2 son dos de las cinco con- 
stantes del movimiento a las que lleva el tratamiento del problema de Kepler 
mediante las variables de accion angulo. 

Necesitamos conocer el promedio de uj en un periodo de la orbita no pertur- 
bada r: 



,., 1 r dAH ^ d fl r ^rr A d{AH) 



Pero el promedio temporal de la Hamiltoniana no perturbada es: 



, . > , , 1 > h f^ dt , ^ 

i^H) = -h{-) = -- -. (26) 



^ /y* 7 L ' /^ / ^ /ya Tc 



Por otro lado, sabemos que / = mr'^(d9/dt), de donde podemos despejar dt 
y sustituirlo en (27), con lo que 
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mh 


(■^^ de 


It . 


/o r"-2 


mh 
It 


\l^ ) Jo 



{AH) = -^ 

r2TT 



ri)]'^-'de . (27) 



donde r] es una fase constante, e es la excentricidad, y donde se ha expresado 
r en funcion de 9 haciendo uso de la ecuacion general de la orbita (con el 
origen en un foco de la conica correspondiente): 



- = -^[l+ecos{e-r])] (28) 



r P 



En el caso en que n = 2: 



(AH) 



2TTmh 



It ' 
(^) = ^. (29) 



En el caso en que n = 3: 



{AH) 



2-^111? hk 



Pt ' 
, , QTTm?hk 



Este ultimo caso, n = 3, reviste una iniportancia especial, ya que la teoria de 
la Relatividad General predice una correccion del movimiento newtoniano 
del orden de r~^ precisamente. Tal prediccion se sometio a prueba con el 
celebre problema de la precesion de la orbita de Mercurio. Sustituyendo 
los apropiados valores de periodo, masa, semieje mayor de la orbita (que 
va incluido en /i), etc., la ecuacion (30) predice una velocidad media de 
precesion: 

(cj) = 42.98 arcsegundos/siglo . 

El valor medido es mucho mayor que el mencionado arriba (por un factor 
mayor que 100). Pero antes de hacer cualquier comparacion deben elimi- 
narse del valor medido, las contribuciones debidas a: a) el efecto conocido 
como la precesion de los equinoccios (movimiento del punto de referenda 
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de longitudes respecto a la galaxia), b) las perturbaciones de la orbita de 
Mercuric debido a la interaccion con los otros planetas. Una vez eliminados 
estos efectos (de los cuales el primero es el de mayor peso), se debe obtener 
lo que seria la contribucion al valor medido de uj, debido al efecto relativista. 
En 1973 se calculo esta ultima contribucion en (41.4ib0.9) arcsegundos/siglo, 
que es consistente con la prediccion que se obtiene de (30). 

9.2 Teoria de perturbaciones independiente del tiempo 

Mientras que la teoria de perturbaciones dependiente del tiempo busca la 
dependecia temporal de los, en un principio constantes, parametros del sis- 
tema no perturbado, la teoria independiente del tiempo pretende hallar las 
cantidades que son constantes en el sistema pertubado. Esta teoria se aplica 
solo a sistemas conservativos y periodicos (tanto en el estado perturbado, 
como en el no perturbado). Por ejemplo, se aplica en el caso de movimien- 
tos planetarios cuando se introduce cualquier perturbacion conservativa al 
problema de Kepler (metodo de von Zeipel o de Poincare). 

Consideraremos aqui el caso de sistemas de un solo grado de libertad. Con- 
sideremos un sistema periodico con una Hamiltoniana indepentiente del 
tiempo de la forma: 



H = H{p,q,X), (31) 



donde A es una constante que especifica la magnitud de la perturbacion y 
se supone suficientemente pequeiia. Asumimos que 



Hoip,q)=H{p,q,0) (32) 



corresponde a un sistema que puede resolverse exactamente (sistema no 
perturbado), por medio del uso de las convenientes variables accion-angulo 
{Jo, Wo): 



Ho{p,q) = Ko{Jo 

I 

'dJ( 



vq = wo = -wy- ; (^0 = t'ot + <5o)- (33) 



La transformacion canonica que nos lleva de {p,q) a {Jo,wo) es independi- 
ente de la forma particular de la Hamiltoniana. Entonces, la Hamiltoniana 
perturbada H{p,q,X) puede escribirse como H{Jo,wo,X). Debido a que la 
Hamiltoniana perturbada si depende de wq, Jq ya no es constante. Por otro 
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lado, en principio, uno puede obtener nuevas variables accion-angulo (J, w) 
apropiadas para el sistema perturbado, tales que: 

H{p,q,X) = E{J,\) 

dE 

J = — — — = 0; [J = constante) . 
ow 

Puesto que la transforniacion que conecta (p, q) a ( Jq, wq) es conocida, ahora 
debemos encontrar la transforniacion canonica S, que conecta {Jq,Wo) a 
{J,w). Si asumimos que A es pequeiia, la transforniacion que buscamos 
no debe diferir niuclio de la transforniacion identidad. La expansion de la 
funcion generadora que buscamos sera entonces: 

S = S{wo, J, A) = So{wo, J) + \Si{wo, J) + \'^S2{wo, J) + ... . (35) 

Para A = requerimos que S sea la identidad; entonces liacemos: 

So = wqJ . (36) 

Las transformaciones canonicas generadas por S son: 

dS . dSi , ^, ,2 ^'S'2 / T^ 

w = -Qj =wo + \-Kj{wo,J) + X -^{wo,J) + ... . 

Jo = ^ = J + X^iwo,J) + X'^iwo,J) + .... (37) 

OWq OWq OWq 

Debido a que wq es una variable angulo (del sistema no perturbado), sabemos 
que At^o = 1 sobre un ciclo. Por otro lado, sabemos que las transformaciones 
canonicas tienen la propiedad de conservar el volumen en el espacio de fase. 
Por tanto, podemos escribir: 



J = <h pdq = <h Jodwo . (38) 

Integrando la segunda ecuacion en (37) sobre una orbita del sistema pertur- 
bado, tenemos: 



(p Jodwo = f Jdwo + ^ A" ® —-^dwo, (39) 
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y sustituyendo (39) en (38): 

dS, 



J = JI^WQ + yX' l^dwo. (40) 

En vista de que Awq = 1) debe cumplirse que: 



n=l 



o bien que: 



/i;^-=°- <^^' 



Ademas, la Haniiltoniana puede expanderse en A como funcion de wq y Jq: 
H{wo, Jo, A) = KoiJo) + XKiiwo, Jq) + X^K2{wo, Jo) + ... , (43) 



donde las Ki son conocidas, ya que H es una funcion conocida de wq y Jq 
para una A dada. Por otro lado tenenios que: 



H{PiqA) = H{wo,Jo,>^) 

= E{J,X), (44) 

que es la expresion para la energia in las nuevas coordenadas de accion 
angulo (donde J sera constante y w sera funcion lineal del tiempo). 

Podemos tambien expandir E en potencias de A: 

EiJ,X) = Eo{J) + \Ei{J) + \^E2{J) + ... . (45) 



En virtud de (44) podemos igualar los coeficientes de las distintas potencias 
de A en (43) y (45). Sin embargo, estas dos expresiones para la energia 
dependen de dos diferentes conjuntos de variables. Para resolver esto ex- 
presaremos Hq en terminos de J, y esto se logra haciendo una expansion de 
Taylor de H{wo, Jq, A) respecto de Jq alrededor de J: 
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H{wo, Jo, A) = H{wo, J, A) + (Jo 



' dj 2 aj2 



+ .. 



(46) 



Las derivadas del desarrollo de Taylor son, hablando con propiedad, derivadas 
respecto a Jq calculadas en Jq = J, si bien podenios escribirlas sin perdidas 
de rigor, en la forma de derivadas con respecto a J, una vez sustituida Jq por 
J en Hq{Jq). En la ecuacion anterior, todo termino que contiene a Jq debe 
reescribirse en terminos de J haciendo uso de la transformacion definida por 
(37) que conecta las coordenadas (Jq, wq) con (J, w). Asi que, de la segunda 
ecuacion en (37) obtenemos (Jq — J), que sustituido en (46) da: 



dH f BSi 



H{wo, Jo, A) = H{wo, «^> A) + — [\^ + A^ 



— + -) + ol77T^ (tt:-) +0(A=^)-(47) 



dw, 



2 9J2 \dwo 



Luego, podemos hacer uso de (43) para expresar H{wo, J, A) = H{wo, Jo, A) | jo=j, 
y sustituir luego en la ecuacion (47) para obtener: 



Hiwo,Jo,X) = Ko{J) + XKiiwo,J) + X^K2iwo,J) + ... 



dwo 



dJ 



dJ 



+ A 



dKo{J)dS2 Id'^Ko 



dSi 



dJ dwo ' 2 9J2 \dwo 
E{J, A) 
^o(^) + A^i(J) + A2^2(J) + .... 



+ ... 



(48) 



Ahora, podemos resolver para los coeficientes Ei{J); esto nos dara la posibil- 
idad de calcular la frecuencia del movimiento perturbado a distintos ordenes 
de perturbacion. Puesto que la expansion de los terminos de Ei no involucra 
una dependencia en wq, entonces la aparicion de wq en (48) debe ser espuria. 
Las Ki{wo, J) de (48) son funciones conocidas, mientras que las Si{wo, J) y 
las Ei{J) son desconocidas. 

Igualando potencias de A tenemos: 



^o(^) 
E,{J) 



+ 



Ko{J) 

Ki{wo,J) + 

K2{wo,J) + 



dSi dKojJ) 
dwQ dJ 
dSidKi{wo,J) 
dwo dJ 



I /dSiY d'^KoiJ) dS2 dKo{J) 



2 \dwr 



5J2 



+ 



dwo dJ 



(49) 
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Venios que, para deterininar Ei necesitamos conocer Si adenias de Ki. 
No debemos perder de vista que la Ei son const antes, pues son solo fun- 
ciones de J (que es una constante del movimiento) . Tambien debemos no- 
tar que dK^/dJ es un termino independiente de wq (pues Kq = Ko(J) = 
Ko{Jo) \jo=j)- Si promediamos sobre wq en ambos lados de la segunda 
ecuacion en (49), obtenemos: 



El = {El) 

Pero hemos visto ya que {dSi/dwo) = §{dSi/dwo)dwo = 0. Luego, 

El = {El) = {Ki) . (51) 

Sustituyendo (51) en el miembro izquierdo de la segunda ecuacion de (49) y 
despejando luego (dSi/dwo), tenemos que: 

dSi^ ^ {Ki) - Ki ^ 

en donde hemos usado que uq = OKq/OJ. 

La solucion para Si se encuentra por integracion directa. En general, vemos 
que el procedimiento para hallar En es el que sigue (una vez asumido que 
se ha resuelto ya para £'„_i): 

• Promediar en ambos lados de la n-esima ecuacion de (49). 

• Insert ar el valor promediado de En que se hallo, en la ecuacion com- 
pleta que se daba en (49) para £'„ (es decir, la que habia antes de 
promediar). 

• La unica expresion desconocida que queda entonces, sera Sn que podra 
obtenerse integrando en: 

dS 

—-^ = funcion conocida de tt;o y ^ . 

• Sustituir Sn en la ecuacion completa para £"„. 
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Una vez hecho esto, se puede continuar para el orden de perturbacion n + 1. 

Segiin podemos ver la determinacion de la energia a un orden particular n 
se determina solo cuando se ha hallado 5„_i, y Sn puede determinarse solo 
cuando se ha hallado En- 

La teoria de perturbaciones independiente del tiempo, es niuy similar al 
esquema de perturbaciones de Rayleigh y Schrodinger en la mecanica ondu- 
latoria. En la teoria ondulatoria, se conoce En solo si la funcion de onda 
se conoce a un orden n — 1. Y ademas, la funcion de onda de orden n se 
encuentra solo cuando se ha calculado En- 
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10. INVARIANTES ADIABATICOS 

Prologo: Un invariante adiabatico es una funcion de los parametros y de las 
constantes del moviniiento de un sistema, que permanece casi constante en el 
limite en el que los parametros cambian infinitamente despacio en el tiempo, 
aunque ellos puedan en ultima instancia cambiar por grandes cantidades. 

CONTENIDO: 

10.1 BREVE HISTORIA 
10.1 GENERALIDADES 
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10.1 BREVE HISTORIA 

La nocion de invarianza adiabatica se remonta a los primeros afios de la 
teoria cuantica. Alrededor de 1910, quienes estudiaban la emision y ab- 
sorcion de radiacion, notaron que los atomos podian existir en estados esta- 
bles en los que su energia estaba fija. Einstein llanio la atencion acerca de la 
cantidad que permanecia casi invariante en un pendulo cuya longitud vari- 
aba continua y lentamente (la cantidad E/u). El sugirio que podrian existir 
cantidades similares asociadas con los sistemas atomicos, que determinaran 
la estabilidad de los niismos cuando estas cantidades adquirieran ciertos 
valores. Tales invariantes adiabaticas fueron luego encontradas por Ehren- 
fest y su uso condujo a la primitiva teoria cuantica de Bohr y Sommerfeld. 
Esta teoria trabajaba bien para los estados del liidrogeno, pero fracasaba 
cuando se le aplicaba a otros atomos. Mas tarde, en 1925-6, surgio una 
teoria cuantica altamente exitosa, debida a Schrodinger, Heisenberg, Born 
y otros; esta teoria empleaba un enfoque diferente. 

El tema de la invarianza adiabatica surgio nuevamente despues de varias 
decadas, en el estudio de iones y electrones en movimiento en el espacio. Este 
tema era del interes de cientificos escandinavos que estudiaban el fenomeno 
de las auroras. Uno de ellos, H. Alfven, mostro en su libro Electrodindmica 
Cosmica que bajo condiciones apropiadas cierta combinacion matematica 
de propiedades de los iones y electrones permanecian constantes a primer 
orden. Aparentemente, Alfven no se percato de que estaba tratando con la 
misma clase de invariantes adiabaticas definidas por Ehrenfest. Fueron L. 
Landau y E. Lifshitz quienes seiialaron la relacion que habia con el tema de 
invariantes adiabaticas. 

10.2 GENERALIDADES 

Mostraremos primero, de una manera sencilla, cual es la cantidad invariante 
adiabatica en el caso del movimiento armonico simple. El metodo a seguir 
consistira en mostrar que, en el limite de la variacion infinitamente lenta 
de los parametros, el invariante adiabatica del sistema unidimensional se 
aproxima a una cantidad que se conserva exactamente en el correspondiente 
sistema bidimensional. 

Consideremos un sistema de un grado de libertad, inicialmente conservativo 
y periodico, que contenga un parametro a que sera inicialmente constante. 
La variacion lenta del parametro no alterara la naturaleza periodica del 
movimiento. Por una variacion lenta entendemos aquella en la que a varia 
ligeramente durante un periodo r del movimiento: 



T{da/dt) <^ a . 
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Pero, aiin cuando las variaciones de a sean pequenas en un periodo cualquiera, 
al cabo de un tiempo suficientemente largo las propiedades del movimiento 
pueden experimental cambios grandes. 

Cuando el parametro a sea constante, el sistema vendra descrito por vari- 
ables accion angulo {wq^Jq) tales que la Hamiltoniana sera H = H{jQ,a). 
Supongamos que la funcion generadora de la transformacion {q,p) -^ (wq, Jq), 
es de la forma: W*{q,'Wo,a). 

Cuando se deja variar a con el tiempo, {wq, Jq) seran aiin variables canonicas 
validas pero W* sera funcion del tiempo a traves de a. Entonces Jq no sera 
constante y wq ya no sera mas funcion lineal del tiempo. La Hamiltoniana 
apropiada sera ahora: 



dW 

K{wQ,jQ,a) = H{Jo,a) + 



dt 
H[Jo,a) + a 



da 



El segundo miembro de la Hamiltoniana se puede ver como una perturbacion 
y la dependencia temporal de Jq viene de: 

. _ dK _ . d ldW*\ 
dwQ dwQ \ da J ' 

Procediendo en una forma analoga a como se procede en la teoria de per- 
turbaciones dependiente del tiempo, buscamos la variacion a primer orden 
en el valor medio de Jq a lo largo del peri odo del movimiento no pertur- 
bado. Como a varia lentamente, podemos considerarla constante durante 
este intervalo y entonces podemos escribir: 



,-, I f . d fdW*\ , a f d fdW*\ , ^,.0 .., 

Puede demostrarse que W* es una funcion periodica en wq, y por consigu- 
iente, tanto ella como su derivada respecto a a, pueden escribirse como una 
serie de Fourier: 



dW* 
da 

Sustituyendo (81) en (80): 



^^^g2^ifc«,o, (5) 
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(Jo) = -- lY,2mkAke^^'^'^Ht + 0{a^,a). (6) 



Como el integrando no tiene ningiin termino constante, la integral se anula. 
Luego, 



(Jo) = + O(d2,a). (7) 



Por tanto, (Jq) no tendra variacion secular de primer orden (o sea en d), 
que es una propiedad deseada de la invarianza adiabatica. De modo que el 
termino casi constante en nuestra definicion de invariante adiabatica, debe 
ser interpretada como constante a primer orden. 
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11. MECANICA DE SISTEMAS CONTINUOS 



Prologo: Todas las forrtiulaciones de la mecanica tratadas hasta ahora han 
estado dirigidas al tratamiento de sistemas que tengan un niimero de grados 
de libertad finito o, como maximo, numerablemente infinito . Sin embargo, 
existen ciertos problemas mecanicos que entranan sistemas continues, como 
por ejemplo el problema de un solido elastico en vibracion. En el cada punto 
del solido continuo participa en las oscilaciones y el movimiento total solo 
puede describirse especificando coordenadas de posicion de todos los puntos. 
No resulta dificil modificar las formulaciones anteriores de la mecanica para 
poder tratar dichos problemas. El metodo mas directo consiste en aproximar 
el sistema continuo a uno que contenga particulas discretas y luego exami- 
nar como cambian las ecuaciones que describen el movimiento, cuando nos 
aproximamos al li'mite continuo. 



CONTENIDO: 

11.1 Formulacion Lagrangiana: de discreto a continuo. 

11.2 Formulacion Lagrangiana para sistemas continuos. 

11.3 Formulacion Hamiltoniana, parentesis de Poisson. 

11.4 Teorema de Noether. 
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11.1 Formulacion Lagrangiana: Transicion de un 
sistema discreto a un sistema continuo. 



Como uno de los casos mas sencillos en el que se puede pasar de un sis- 
tema discreto a uno continuo consideremos el de una varilla elastica in- 
finitamente larga que efectiia pequefias vibraciones longitudinales, es decir, 
desplazamientos oscilantes de las particulas de la varilla paralelos a su eje. 
Un sistema compuesto por particulas discretas que se aproxime a la varilla 
continua es una cadena infinita de puntos materiales iguales separados por 
distancias a y unidos por resortes uniformes sin masa de constante de rigides 
k (ver la figura). 



despiazado 
del 



Supondremos que los puntos materiales solo pueden moverse a lo largo de la 
direccion de la cadena . Podemos ver que el sistema discreto es una extension 
de la molecula poliatomica lineal tratada en el capitulo 6 de Goldstein. Po- 
dremos, pues, obtener las ecuaciones que describen el movimiento mediante 
las tecnicas habituales para oscilaciones pequefias. Representando por r/, el 
desplazamiento de la particula i-esima respecto a su posicion de equilibrio, 
la energia cinetica es 



t=Iy.^vI 



:i) 



donde m es la masa de cada particula. La energia potencial correspondiente 
es la suma de las energi'as potenciales de cada resorte a consecuencia de 
hallarse estirado o comprimido respecto a su longitud natural: 



1 



^ = 7j zi ^ (^i+1 ~ ^^y 



(2) 



De las ecuaciones (|l|) y (Q) obtenemos la Lagrangiana del sistema 



1 ^ / .2 , . ,2 



L = T-V=-J2(mV, -A:(7?,+i-r/,) 



(3) 



que tambien puede escribirse de la forma 



^mz' 



- '1^ -kg "•+■ - "' 
a \ a 



^aLi, 



(4) 
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donde a es la separacion de equilibrio entre puntos. Las ecuaciones de 
movimiento de Lagrange para las coordenadas r/j resultan ser 

La forma particular de L de la ecuacion (Q) y las ecuaciones de movimiento 
correspondientes las hemos elegido por ser comvenientes para efectuar el 
paso al limite a una varilla continua al tender a cero a. Esta claro que m/a 
se reduce a la masa por unidad de longitud ^ del sistema continuo, pero el 
valor limite de ka no resulta tan evidente. Recordemos que en el caso de 
una varilla elastica que cumpla la ley de Hooke,el alrgamiento de la varilla 
por unidad de longitud es directamente proporcional a la fuerza o tension 
ejercida sobre ella, relacion que podemos escribir en la forma 

donde ^ es el alargamiento por unidad de longitud e y es el modulo de Young. 
Ahora bien, el alargamiento de una longitud a de un sistema discreto, por 
unidad de longitud, sera ^ = (ryi+i — rji) /a. La fuerza necesaria para estirar 
el resorte esta cantidad es 



k (?7i+i - ?7i) = ka{ -^ j 



por lo que Ka debe corresponder al modulo de Young de la varilla continua. 
Al pasar del caso discreto al continuo, el indice entero i que identifica el punto 
material particular se convierte en la coordenada de posicion continua x; en 
vez de la variable rji tenemos rj (x) . Ademas, la cantidad 

m+i - m _ vix + a) -r]{x) 
a a 

que figura en L, tiende evidentemente al limite 

dr] 
dx' 

cuando a tiende a cero. Por ultimo, la suma extendida a un mimero discreto 
de particulas se convierte en una integral extendida a x, la longitud de una 
varilla, y la Lagrangiana (H) queda en la forma 



i/(-'- 



L = ^I{M -Y[^) ]dx. (6) 



En el limite, cuando a tiende a cero, los dos liltimos terminos de la ecuacion 
de movimiento (^) resultan ser 



Lima, 




dr] 
dx 
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tomando de nuevo el limite cuando a tiende a cero la ecuacion define clara- 
mente la segunda derivada de r]. Por tanto, la ecuacion de niovimiento para 
la varilla elastica sera 

que es la conocida ecuacion de onda en una dimension con velocidad de 
propagacion 

^ (8) 

La ecuacion (P) es la conocida formula de la velocidad de propagacion de 
las ondas elasticas longitudes. 

Este sencillo ejemplo es suficiente para ilustrar las caracteristicas principales 
de la transicion de un sistema discreto a uno continuo. El hecho mas impor- 
tante que hemos de comprender es el papel que desemperia la coordenada 
X. No se trata de una coordenada generalizada; solo hace las veces de indice 
continuo que sustituye al indice discreto i. Al igual que cada valor de x 
corresponde una coordenada generalizada rj (x) . Como ij depende tambien 
de la variable continua t, debemos tal vez escribir con mayor precision que 
r] (x.t). Indicando que x, al igual que t, puede considerarse como parametro 
que entra en la lagrangiana. Si el sistema continuo fuese tridimensional y 
no unidimensional, como en este caso, las coordenadas generalizadas se dis- 
tinguirian mediante tres indices continuos x,y,z y se escribirian en la forma 
r] {x,y,z,t) . Notemos que las cantidades x,y,z,t son totalmente indepen- 
dientes unas de otras y solo aparecen en t] como variables explicitas. Las 
derivadas de r] respecto a cualquiera de ellas podran, pues, escribirse siem- 
pre en forma de derivadas totales sin ninguna ambigiiedad. La ecuacion 
@ indica tambien que la lagrangiana aparece como integral para el indice 
continuo x; en el caso tridimensional, la lagrangiana tendria la forma 

L= f f f Cdxdydz, (9) 

donde C se denomina densidad Lagrangiana. En el caso de vibraciones 
longitudinales de la varilla continua, la densidad Lagrangiana es 

2 / ^^ \ 2 ~i 



drA _y(^\ 



y corresponde al limite continuo de la cantidad Lj que aparece en la ecuacion 
@. Es la densidad de Lagrangiana, mas que la propia Lagrangiana, la que 
utilizaremos para describir el movimiento del sistema. 



11.2 Formulacion Lagrangiana para sistemas con- 
tinuos 

Notemos en la ec. @ que la C para la varilla elastica depende de V= drj/dt, 
la derivada espacial de rj, drj/dx; x y t desempefian un papel similar al 
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de los parametros de esta. Si ademas de las interacciones entre vecinos 
mas proximos hubiesen fuerzas locales, L fuera funcion de la misma t]. En 
general £ para todo sistema continue, puede ser funcion explicita de x y 
t. Por tanto, la densidad de Lagrangiana para todo sistema continuo debe 
aparecer para todo sistema continuo unidimensional de la forma 

La Lagrangiana total siguiendo la forma de la ec.(|9|) sera 

L = Cdx, 
y el principio de Hamilton en el limite del sistema continuo adopta la forma 

61 = 6 f f Cdxdt = 0. (12) 



1 

Del principio de Hamilton para el sistema continuo, debera ser posible de- 
ducir el limite continuo de las ecuaciones de movimiento, para esto como en 
la seccion 2-2 de Goldstein podemos obtener en el espacio rj un camino vari- 
ado de integracion conveniente, eligiendo rj entre una familia de funciones 
de T] dependiente de un parametro: 

rj {x,t;a) = r]{x,t;0) + aC {x,t) . (13) 

Donde r] {x,t;0) esla funcion correcta que satisface el principio de Hamilton 
y C es una funcion cualquiera de buen comportamiento que se anule en los 
puntos extremos en t y en x. Si consider amos I funcion de a, para que sea 
una extremal para la derivada de / respecto de a se anulara en a = 0. Ahora 
por derivacion directa de I tenemos 



/"*' r'w rl /^^^'? ^^ ^ f^^\ ( dC\ d (dr]\ 



dl_ 

da Jti Jxi """""' I 07? 9a ' d^da\dt, ^^ ^ 

(14) 
Como la variacion de rj, es decir a^, se anula en los puntos extremos, inte- 
grando por partes segiin x y t obtenemos las relaciones 

*2 dC d fd7]\ [^2 d ( dC 



y 

'■^a dC d /dr]\ , /"^a d 

dx 





ii d^da\dxj Jxi dx 

De aqui el principio de Hamilton se podra escribir de la forma 



t2 
tl 



^2 , , fa£ d ( dC\ d ( dC\\ / di]\ 
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Ahora debida a la naturaleza arbitraria del caniino variado implica que la 
expresion entre Haves es cero: 

La ecuacion anterior corresponde a la ecuacion correcta de movimiento de- 
ducida del principio de Hamilton. 

En el caso concrete de las vibraciones longitudinales en una varilla elastica, 
la forma de la densidad de Lagrangiana ec. (|lO| ) indica que 

dC dr] d£ drj dC 

'Q^~^~di' 'd^~~ ^' '&n ~ ' 

Asi pues, tal como queri'amos, la ecuacion de Euler-Lagrange (|^) se reduce 
a la ecuacion de movimiento (j^)- 

La formulacion de Lagrange que acabamos de desarrollar corresponde a sis- 
temas continues, evidentemente, se pude generalizar a sistemas bi- , tridi- 
mensionales y cuadridimensionales. Matematicamente, conviene pensar en 
un espacio cuadridimensional de coordenadas Xq = t,xi = x, X2 = y,xs = z. 
Para un manejo matematico mas facil introducimos la siguiente notacion 

^ dij^ ^ drj_ ^ d'^rg 

Con esta notacion, y haciendo la extension a un espacio cuadridimensional, 
la forma general de la densidad de Lagrangiana (|T^) toma la forma: 

£ = £(r/p,r/p,^,§^). (18) 

La Lagrangiana total es entonces una integral extendida al espacio tridimen- 
sional: 

L= f C{dxi). (19) 



Al principio de Hamilton corresponde una integral extendida a una region 
de un espacio cuadridimencional 

61 = 6 f C{dx^)=0, (20) 

donde la variacion de las ijp se anulan en la superficie S de contorno de 
la region de integracion. La deduccion de las correspondientes ecuaciones 
de movimiento de Euler-Lagrange tiene lugar simbolicamente como antes. 
Consideremos un sistema de funciones 

Tjp {xy\ a) = rjp (xu) + aC (x^) 
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variadas de un solo parametro que se reducen a rip [xu) cuando el parametro 
a tiende a cero. La variacion de / equivale a hacer igual a cero la derivada 
de I respecto a a es decir: 



dl_ 

da 



dC dr]p dC drjp^^ 



di]p da drjp^y da 
Integrando por partes la ec. (^), tenemos 



{dxp 



0. 



(211 



dl 
da 



dC 

drip 



d 



dC 



dxy \dripy 



drip 
da 



idxp)+ / (dxp) 



d 



dC drip^i, 
dVp,u da 



y tomando el limite cuando a tiende a cero la expresion anterior se reduce 
a: 



dl\ 

da J I 



{dx^ 



dC 

drip 



d 

dXy 



dc 

dn 



'p,v 



KdajQ 



(22) 



Ahora debido a la naturaleza arbitraria de la variacion de cada tipo de rip 
significa que la ecuacion (g^) cuando sea nulo cada uno de los corchetes por 
separado esto es: 

d ( dC \ _d£_ 
dxy \di!ip^y) drip 

Las ecuaciones (53) representan un sistema de ecuaciones en derivadas par- 
ciales para las cantidades campo, que tiene tantas ecuaciones cuantos valores 
diferentes de p haya. 

Ejemplo: Dada la densidad Lagrangiana para un campo acustico 



£ = - ( ^0 r/ +2PoV • rf- 7P0 (V • 7?)' 



/io es la densidad masica de equilibrio y Pq la presion de equilibrio del gas. 
El primer termino de C es la densidad de energia cinetica, mientras que los 
terminos restantes representan el cambio que sufre la energia potencial del 
gas por unidad de volumen a consecuencia del trabajo efectuado sobre el 
gas o por el curso de las contracciones y expansiones que son la marca de 
las vibraciones aciisticas, 7 es el cociente entre los calores molares a presion 
y a volumen constante obtener las ecuaciones de movimiento. 

Solucion: 

Con la notacion cuadridimensional, la densidad de Lagrangiana queda en la 
forma 

1 



^ = 2 if^OVi,OVi,0 + '^PoVi,: 



lPoVi,ir]j,j) ■ (24) 

De la ecuacion ( [23| ) se obtienen las siguientes ecuaciones de movimiento 

/^o??i,oo - iPoVuj = 0, j = 1, 2, 3. (25) 
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Volviendo a la notacion vectorial, las ecuaciones (25) toman la forma 

(f 
Ho—dt^ - 7P0VV ■ff=0. (26) 

rj 

Ahora utilizando el hecho de que en vibraciones de pequena amplitud la 
variacion relativa de la densidad del gas esta dado por la ecuacion 

a = —V • 7/ . 

Ahora aplicando la divergencia y utilizando la ecuacion anterior obtenemos 

7P0 dt^ 
la cual es una ecuacion de onda tridimensional, siendo 



V /^o 
la velocidad del sonido en los gases. 

11.3 Formulacion Hamiltoniana, parentesis de Pois- 
son. 

11.3.1 Formulacion Hamiltoniana 

La formulacion de Hamilton para sistemas continuos se hace en forma pare- 
cida como se hace para sistemas discretes. Para indicar el procedimiento 
volvamos a la cadena de puntos materiales tratada anteriormente, donde 
para cada r]i hay una cantidad de movimiento canonica 

dL dLi 

Pi = — — = a — -. (27) 

dVi dr], 

La Hamiltoniana del sistema sera, entonces 

H = Pi Vi -L = a—^ Vi -L, (28) 

dVi 

o sea 

H = a(^%-L,) . (29) 

\dVi J 

Recordando que en limite cuando a tiende a cero, L — > £ y la suma de la 
ecuacion (|2^) se convierte en una integral por lo que el Hamiltoniano toma 
la forma: 

^ = /'^-(^^-^) • (30) 
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Las cantidades de movimiento canonicas individuales pi, dadas por la ecuacion 
(p7|), se anulan en el limite de la continuidad, pero podemos definir una den- 
sidad de cantidad de movimiento vr que permanezca finita: 

Lim^^o— = TT = --^ . (31) 

a dV 

La ecuacion pO| ) tiene la forma de integral espacial de una densidad de 
Hamiltoniana 7i definida por 

n = 7rV -C . (32) 

Aun cuando se pueda introducir asi una formulacion de Hamilton de man- 
era directa para campos clasicos, demonos cuenta de que el procedimiento 
singulariza la variable tiempo a la que habra que darle un tratamiento espe- 
cial. Contrasta, pues con el desarrollo que hemos dado a la formulacion de 
Lagrange en el cual se trataban simetricamente las variables independientes 
del tiempo y espaciales. Por esta razon, el metodo de Hamilton, se tratara 
en forma un tanto distinta. 

La via evidente para la generalizacion a un campo tridimensional descrito 
por cantidades campo r]p es la siguiente: 
Se define una cantidad de movimiento canonica 

-. K) = 1^- (33) 

Donde las cantidades r]p {xi,t) ,-Kp {xi,t) juntas, definen el espacio fasico de 
infinitas dimensiones que describe el campo clasico y su desarrollo en el 
tiempo. 

Analogamente como en el sistema discrete podemos hallar un teorema de 
conservacion para vr que sea algo parecido al correspondiente al de cantidad 
de movimiento canonico de sistemas discretos. Si una cantidad de campo 
dada r]p es ciclica es decir que C no contenga explicitamente a r/p, la ecuacion 
de campo de Lagrange presenta el aspecto de enunciado de la existencia de 
una corriente conservativa: 

d dC 



o sea 







^-^4^ = 0. (34) 

dt dx. Of] ■ 

Se sigue que si es ciclica ijp, existe una cantidad integral que se conservativa 

Up = J dVTT^{Xi,t). 

La generalizacion para la densidad de la ecuacion (32) para una densidad 
Hamiltoniana es 

[Vp, r]p,, , T^p: Xf}j = TTp Vp -C, (35) 



n(j 
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donde se supone que se puede eliminar la dependencia funcional de Vp por 
inversion de las ecuaciones de definicion (p3). De esta definicion se deduce 
que 



on 



--V„ +vr. 



dv. 



dC dV. 



dir^ drj^ Ott^ 



--v„ 



(36) 



en virtud de la ecuacion (^3[). Analogamente obtenemos 



dn 

drjp 



dr]. 



dC dV, 



TT, 



dri d'h. di] 



dC 
dri 



dC 



(37) 



Ahora utilizando las ecuaciones de Lagrange en la ec.(37) queda como 

' d£ \ . d ( dC\ 



dn 

di]p 



_d_ 



dr]^ 



vr 



d 
dx. 



.^^P. 



(38) 



Debido a la aparicion de C aun no tenemos una forma util. Sin embargo, 



haciendo una deduccion analoga a la de los terminos j^ — y "" — 



tenemos 



dn 

dn 



dr]^ 



dC dV, 



dC 



vr. 



'p,i 



di]^^ d V^ dr]^^ drj^^ 



dC 



Por lo tanto sustituyendo ( p9| ) en (38) obtenemos 



dn 

drjp 



d 
dxi 



dn 

drjp^i 



IT 



(39) 



(40) 



Las ecuaciones (36) y (4C) las podemos expresar con una notacion mas 
proxima a la de las ecuaciones de Hamilton para un sistema discrete intro- 
duciendo la nocion de derivada funcional definida en la forma 



d 



d d 



5tp dtp dxi dip^i 



Como n no es funcion de vr^ . las ecuaciones ( |3' 



en la forma 



5n 

6tt„ ' 



vr 






(41) 
^ se pueden escribir 

(42) 



Ahora con la misma notacion las ecuaciones de Lagrange (p3| ) toman la 
forma 

^f^^-f.O, (43) 

Sin embargo, la ventaja casi linica de la derivada funcional estriba en la 
semejanza resultante con un sistema discreto. Por otra parte, sorprende el 
tratamiento paralelo de las variables temporal y espaciales. 
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11.3.2 Parentesis de Poisson 

Podemos obtener otras propiedades de TL desarrollando la derivada total 
respecto al tiempo de la ecuacion (^5|), recordando que hay que considerar 
que ?7p es funcion de ilp^Vpjj'^p Y ^m- Tenemos, entonces que 



dn 


dv 


dc . 


dC dr] 


dC dr]. 


dc 




=7r 1] +7r ^ - 


v 






— 


dt 


" " " dt 


dvp " 


dV^ dt 


dVp, dt 


dt 



En la expresion el segundo termino y el cuarto se aniquilan debido a la 
definicion (|33|), por lo que la derivada se simplifica quedando 

dn . ■ dC . dC drj^^ dC 

=7r Tj r] '-PiL . (44) 

dt " " dr]p d7]^^ dt dt ^ ' 

Por otra parte, considerando TC funcion de 77^, r/ . , vr^ y vr^, la derivada total 
respecto al tiempo es 

(£{_. dn dn ■ dn_dr]^ dn 
~dF~'''d^p^d^p'^'^dv^i^^~dt' ^ ' 

donde la expresion se escribio de tal manera que facilite la comparacion con 
el segundo miembro de la ecuacion (0), donde usando las ecuaciones (|3^), 



(37) y (pgp obtenemos 

-dT = -W ^''^ 

la cual es analoga a la correspondiente para sistemas discretos. 
En canibio no se cumple que las derivadas total y parcial respecto al tiempo 
no son, en general iguales. Utilizando las ecuaciones de movimiento de 
hamilton (ec. (^) y (^0[)) e intercambiando los ordenes de derivacion, la 
ecuacion (|4^) se puede escribir como 

dn _ dn d f dn\ dn dv^ dn 

dt d-Kp dxi \ di]p^i I drjp^i dxi dt 



Ahora utilizando la ecuacion (|4q ) y combinando terminos tenemos final- 
mente 

dn__d_(. dn_\ dn 
~dF~d^iy''d^i)^~dt' ^ ' 

que es lo mas que podemos aproximarnos a la ecuacion correspondiente para 
sistemas discretos. 

Por otra parte cuando C no contenga a t explicitamente, tampoco la con- 
tendra n esto implica la existencia de una corriente consservativa y por lo 
tanto la conservacion de una cantidad integral, en este caso 

H = fndV . (48) 



130 



Asi pues, si TL no es funcion explicita del tiempo, la cantidad que se conserva 
no es "H, sino la cantidad integral H. 

La Hamiltoniana no es mas que un ejemplo de funciones que son integrales 
de volumen de densidades. Podemos formular directamente un formalismo 
general para la derivada respecto al tiempo de dichas cantidades integrales. 
Consideremos una cierta densidad U que sea funcion de las coordenadas del 
espacio fasico {rjp,TTp), de sus gradientes espaciales y posiblemente de x^: 

U=u{r]p,iT^,r]p^i,TT^^^,Xpj . (49) 

La cantidad integral correspondiente es 

U{t) = fudV (50) 



donde la integral de volumen se extiende a todo el espacio limitado por la 
superficie de contorno sobre la cual se anulan rjp y TTp. Derivando U respecto 
al tiempo tenemos en general, 

dU [ { dU . dU . dU . dU . dU\ 

Consideremos un termino tal como 

dU ■ r „BU dV 



dV^ Vo.i= dV 



p 



dr]p,i ' J dr]p^i dxi 

Ahora intengrando por partes, considerando que r^p y las derivadas se anulan 
en las superficies de contorno, tenemos 

dV-^ ' =-[dV' —(-^\ 

drjp^i '''' J " dxi ydrjp^ij ' 

Para el termino en Tip j se hace un procedimiento similar. Sustituyendo las 
expresiones obtenidas y agrupando coeficientes de ^ y de TTp respectivamente, 
y usando la notacion de derivada funcional la ecuacion (^TI) se reduce a 

dU r ._ ( 6U . 5U . dU' 



Por ultimo, introduciendo las ecuaciones de movimiento canonicas ( |4^ ), 

tGIlGIXlOS 

^=/,y(^i!^-^i?^U/,y^. (53) 

dt J y dr]p diTp orjp diTp j J ot 

La primera integral del segundo miembro corresponde claramente a la forma 
de corchete de Poisson. Si U y W son dos funciones de densidad, estas 
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consideraciones nos sugieren la definicion del corchete de Poisson de las 
cantidades integrales como 

lV,Wl^f^\fi'^-'^fi]. (54) 

J [ drip dvTp dr]p diTp J 

Definamos tambien la derivada parcial de U respecto a t, mediante la sigu- 
iente expresion 

La ecuacion ( |5^ ) podra entonces escribirse en la forma 

que corresponde precisamente, en esta notacion a la ecuacion para sistemas 
discretos. Como por definicion, el corchete de Poisson de H consigo misma 
es nulo, la ecuacion (|4q ) se concretara en 

^ - ^ (57) 

dt - dt ' ^^^> 

que es la forma integral de la ec. (^). Asi pues, el formalismo de corchetes 
de Poisson aparece como consecuencia de la formulacion de Hamilton. Pero 
no podemos llevar a cabo una descripcion por corchetes de Poisson de la 
teoria de campos en correspondencia paso a paso con la de los sistemas 
discretos. 

Sin embargo, hay una manera de tratar los campos clasicos que provee 
casi todo lo de la formulacion de Hamilton y de corchetes de Poisson de la 
Mecanica para sistemas discretos. La idea fundamental de este tratamiento 
es sustituir la variable espacial continua o el indice continuo por un indice 
discrete numerable. 

El requisito de que 77 se anule en los extremos es una condicion de contorno 
que se podria realizar fisicamente colocando la varilla entre dos paredes per- 
fectamente rigidas. Entonces, la amplitud de oscilacion se puede representar 
mediante una serie de Fourier: 

v{x) = Yl 1^ ^^^ — ^ — • (^^) 

n=0 

En vez del indice continuo x tenemos el indice discreto r/. Podremos utilizar 
esta representacion de x solamente cuando r/ [x) sea una funcion regular, 
cosa que sucede en la mayoria de cantidades de campo fisicas. 
Supondremos que solo hay una cantidad campo real r] que se puede desar- 
rollar en serie de Fourier tridimensional de la forma 

^ ^^' *) = T7T72 H 9^= (*) ^^P y^ ■ ^) (^9) 

fc=0 
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Aqui, k es un vector de onda que solo puede tomar modulos y direcciones 
discretos de manera que en una dimension lineal dada solo encaje un numero 
entero (o a veces semientero) de longitudes de onda. Decimos que k tiene un 
espectro discrete. El subindice escalar k representa una cierta ordenacion del 
sistema de indices enteros que se utiliza para enumerar los valores discretos 
de k; V es el volumen del sistema, el cual aparece en forma de factor de 
normalizacion. 

La ortogonalidad de las exponenciales en todo el volumen se puede enunciar 
mediante la relacion 

1 /■ i{k-k 
V 



e ^ ^ dV = 6,y . (60) 



En realidad, los valores permitidos de k son aquellos para los cuales se 
satisface la condicion (pO|), y los coeficientes q^ (t) estan dados por 



qk{t) = ^Je~^^-^7^{r,t)dV. (61) 

De manera analoga para la densidad de cantidad de movimiento canonica 
tenemos 

^{^,t) = ^Y.Pkit)e-'~^-^ (62) 

k 
con pk (t) definido como 

Pk{t) = ^Je-^^-^7,{y,t)dV. (63) 

Tanto Qk como pk son cantidades integrales. Podemos, pues, buscar los 
corchetes de Poisson de dichas cantidades. Como las exponenciales no con- 
tienen las cantidades campo tenemos, por la ecuacion (^) 

[Qk,Pk'] ■ 



^ 1 dVe-^'^-^l^''^^ - 


6tt 6r] 
6r} 6n 


l^JdVe-'^-^ 





o sea, por la ecuacion (pQ), 

[qk,Pk'] =K,k' ■ (64) 

De la definicion de los corchetes de Poisson resulta evidente que 

bk^Qk'] = [Pk,Pk'] =0 • (65) 

La dependencia temporal de qk se hallara a partir de 

qu it) - [qk,H] -^^,JdVe |-— - — - 
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o sea 

(66) 



(67) 
(68) 







9fc (t) = 


yi/2 J ' 


„-i k -T* 


6n 

Sir 


Por otra parte, 


tenemos que 












dH 

dpk 


^JdV 


dn dpk 




por lo que obtenemos 


dir 


1 


-i k -T* 





dpk FV2 

Comparando las ecuaciones ( |67| ) y ( p^ ) tenemos 

■ , , dH 

De manera similar podemos obtener la ecuacion de movimiento para pk 

Pu= -^ . (70) 

dqk 

Asi las cantidades Pk Y Qk: obedecen pues, las ecuaciones de movimiento de 
Hamilton. 



11.4 Teorema de Noether 

Ya hemos visto multiples veces que las propiedades de la Lagrangiana (o de 
la Hamiltoniana) implican la existencia de cantidades conservativas. Asi, 
si la Lagrangiana no contiene explicitamente una coordenada particular 
de desplazamiento, se conserva la correspondiente cantidad de movimiento 
canonica. La ausencia de dependencia explicita de la coordenada significa 
que la Lagrangiana no queda afectada por una transformacion que altere el 
valor de dicha coordenada; se dice que es invariante o es simetrica ante la 
transformacion dada. 

La simetria ante una transformacion de coordenadas se refiere a los efectos 
de una transformacion infinitesimal de la forma 

donde la variacion 6x^ puede ser funcion de las demas x^- El teorema de 
Noether considera tambien el efecto de una transformacion de las propias 
cantidades campo, la cual podemos escribir en la forma 

V (Xf,) -^ Vp (x^,) = Tip {Xp) + Slip {Xp) ■ (72) 

Aqui 5rip (x^) mide el efecto de las variaciones de x^ y de 77^ y puede ser 
funcion de las demas cantidades campo r/^ . La variacion de una de las vari- 
ables campo en un punto particular del espacio x^ es una cantidad diferente 

% (x^ = Vp {xj + '67]^ {x^) . (73) 
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La descripcion de las transformaciones en funcion de varaiaciones infinitesi- 
males a partir de las cantidades no transforniadas nos indica que solo estamos 
tratando con transformaciones continuas. Asi la simetria ante la inversion 
en tres dimenciones no sera una de las simetrias a las que se pueda aplicar 
el teorema de Noether. A consecuencia de las transformaciones tanto de 
las coordenadas como de las cantidades campo, la lagrangiana aparecera, en 
general, como funcion diferente de las variables campo y de las coordenadas 
del espacio y tiempo: 

-C {rip (^ J , %,. {x^) , x^) -^ C' (ji^ (x'^) , rj^^ (x^) , x^) . (74) 

La version del teorema de Noether que vamos a presentar no constituye 
la forma mas general posible pero facilita la deduccion sin restringir de 
manera importante el ambito de aplicacion del teorema ni la utilidad de las 
conclusiones. Supondremos que se cumplen tres condiciones: 

1. El cuadriespacio es euclideo. Este requisito, restringe el espacio-tiempo 
relativista al espacio de Minkowski, que es complejo pero euclideo. 

2. La densidad de lagrangiana presenta la misma forma funcional para las 
cantidades transformadas que para las cantidades originates, es decir, 

3. La magnitud de la integral de accion es invariante ante la transfor- 
macion, es decir, 

l' = J^, {dxj C' (t]'^ (x^) , 77^^ (x^) , x'^) = j^C {t]^ (x J , 7]^^ (xj , x^ 

(76) 
La combinacion de las ecuaciones (|75|) y (^) nos da el requisito 

.. (dXp) C {ri^ (x J , 77^^ (xj , x^) - / £ (r/^ (xj , 77^,, (x J , xj = . 

(77) 
De la condicion de invariancia, la ecuacion (^Tj) adopta la forma 

/ , dx^C (77', x^) - / (ix^£(77,xj (78) 

^(V,a;^) - ^{v.x^) + / C{r])6x^dS^, = . 



dx„ 



n 



Aqui C (77, x^) es una abreviatura de la dependencia funcional total, S es la 
superficie tridimensional de la region 0, y 5x^ es de hecho el vector diferencia 
entre puntos de 5 y los puntos correspondientes de la superficie transfor- 
mada S . La ultima integral se puede transformar mediante el teorema de 
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la divergencia cuadridimensional con lo que tendremos para la condicion de 
invarianza 

= l^dxA[c[r]' ,x^) - C{r^,x^)] + ^ {C{7],x^) 6x^)\ . (79) 

Ahora usando la ecuacion (|73|) , el termino entre los corchetes puede escribirse 
en primera aproximacion en la forma 

^ (i K) ' V. (a;J ,3;^) --C (??, {xj,v,,. (xj ,a;J = —6r,^ + a^pH.^- 
Utilizando las ecuaciones de campo de Lagrange 

CW,xJ-C{r,,x^) = ^^(^-^-6rj^y 
Luego la condicion de invarianza ( [79[ ) aparece en la forma 

(80) 



i^'^'^M^'''-"^']'"' 



que ya tiene la forma de una ecuacion de corriente conservativa. 
Sin embargo resulta litil desarrollar algo mas la condicion especificando la 
forma de la transformacion infinitesimal en funcion de R parametros in- 
finitesimales Err = 1,2,...,R, tales que las variaciones de x^ y r]^ sean 
lineales en los £r'- 

6x^ = ErXru, St]^ = 6^*,,^. (81) 

Sustituyendo estas condiciones en la ecuacion (|80|) obtenemos 



dx^ [\dr]^^ "•" ""J dr]^^ 

Como los parametros £r son arbitrarios, existen r corrientes conservativas 
con teoremas de conservacion diferenciales: 

" ^^^'^-,.„-«„U..-^*.,l=0. (82) 



dxu 1 \dVn,. "'" ""/ dr] 



rp 
P,' 



Las ecuaciones (82) constituyen la principal conclusion del teorema de Noether, 
el cual dice pues, que si el sistema tiene propiedades de simetria tales que se 
cumplan las condiciones (1) y (2) para transformaciones del tipo indicado 
en las ecuaciones (^), existiran r cantidades conservativas. 
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